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EDITORIAL 


A regularidade que se pretendia para a revista foi conseguida. 
Este é o terceiro número desta nova fase e o primeiro, já, de 1991. 


Consideramos extremamente positivo que estes factos se tenham acompanhado de 
um imediato aumento de interesse e de colaboração: originais já enviados, 
promessa de outras valiosas colaborações, artigos convidados que já são uma 
realidade num futuro próximo. Isto mostra que os docentes do IST, e ainda outros, 


estavam atentos e punham expectativa nos rumos da TÉCNICA. 


No entanto, o original de que dispomos, embora valioso e não escasso, ainda não é 


suficiente para garantir um crescimento automático e natural da revista. 


=“ 


Outro factor material de limitação é a publicidade veiculada pela revista. 
Julgamos que no universo vasto que a TÉCNICA interessa, haverá quem encontre 
motivos fortes para anunciar na revista. Além da acção normal de qualquer 
publicidade, representa um investimento importante nos docentes mais jovens e 


ainda nos futuros engenheiros. 

Finalmente, desejamos salientar que sem o apoio importante do Conselho 
Directivo do IST não teria sido possível conseguir o que em tão pouco tempo se 
conseguiu. Apreciávamos, era justo, que outros seguissem o exemplo do Conselho 


Directivo desta Escola. 


2 | Técnica - 91/1 


NORMAS DE ADMISSÃO DE TRABALHOS E 
INSTRUÇÕES PARA OS AUTORES 


1 A revista TÉCNICA, em princípio, 


publicará anualmente quatro números de 
natureza técnico-científica. 


Atendendo ao contexto onde se insere e 
público a que, em conformidade e de modo 
principal, se destina, a revista integra-se a 
partir das três seguintes formas de 
publicação: 


ARTIGOS POR CONVITE 

São da iniciativa da Direcção da Técnica e 
destinam-se, em princípio, a trabalhos de 
síntese no domínio das ciências básicas ou 
ciências da engenharia. Além do mais, 

rocura-se que estes artigos se constituam 
ontes recomendadas de estudo aos alunos 
e de consulta para investigadores e 
docentes. 


ARTIGOS TECNICO-CIENTIFICOS 
São artigos submetidos a publicação. 
Destinam-se essencialmente a dar noticia 
sobre aquisições actuais da ciência e da 
técnica mas em termos que se vá além da 
simples divulgação, sendo desejável que 
nesses trabalhos figurem parâmetros e 
referências que permitam elaborar juízos 
críticos e comparativos. Fundamental- 
mente, uma notícia sucinta, com carácter 
crítico, dando conta de uma aquisição 
recente e importante para a técnica é a 
marca fundamental destes artigos que vão 
traduzir, em princípio, investigação 
científica realizada pelo autor. 


NOTAS CIENTIFICAS 

Destinam-se essencialmente à 
apresentação de novos resultados no 
domínio da investigação científica ou do 
desenvolvimento tecnológico. Pretende-se 
uma nota curta, realçando os resultados 
alcançados e referenciando conveniente- 
mente trabalhos afins, evitando a 
exposição longa de antecedentes, que 
deverá ser remetida para referências. 
Deste modo, a originalidade que se 
pretende e se traduzirá pelos resultados 
novos que o autor apresenta é a marca 
fundamental de uma nota científica. 


er que tal se ofereça conveniente, a 
TÉCNICA ouvirá especialistas sobre a 
essência e estrutura dos trabalhos subme- 
tidos a publicação, com natural excepção 
dos Artigos por Convite. Quando o autor ou 
autores forem alunos é necessário que do 
trabalho conste a indicação de um super- 
visor de reconhecida idoneidade científica e 
técnica que tenha comparticipado na 


análise crítica e prévia discussão do 
trabalho em causa. 


Os Artigos por Convite e os Artigos 
Técnico-Científicos, deverão ser escritos em 
língua portuguesa. As Notas Científicas 
poderão ser escritas em inglês ou francês. 


Os autores deverão atender na elaboração 
dos artigos às seguintes condições: 

- os Artigos por Convite não poderão 
exceder 30 páginas dactilografadas, in- 
cluindo o espaço reservado a ilustrações; 

- os Artigos Técnico-Científicos não 
epi exceder 15 páginas dactilogra- 
adas, incluindo o espaço reservado a 
ilustrações; 

- as Notas Científicas não poderão exceder 
8 páginas dactilografadas, incluindo o 
espaço reservado a ilustrações; 

- nome exacto com que desejam figurar; 

— ordem de entrada dos autores, havendo 
mais do que um; 

—- instituição(ões) a que pertençam e que 
desejem referir, indicando claramente o 
nome dessa instituição, morada e sigla. 


Os trabalhos devem ser dactilografados a 
dois espaços, em folhas modelo A4, 
devidamente numeradas. 


Os trabalhos devem ser enviados em 
triplicado, recomendando-se aos autores 
que conservem uma cópia. 


Conjuntamente com o trabalho, devem ser 
enviados: 
— título do trabalho em português e inglês; 
- resumo em português e abstract em 
inglês, pelo menos; 
- nome, direcção e telefone do autor e 
demais dados que facilitem o processo de 
revisão de provas. 


Os autores deverão assegurar-se, antes de 
enviarem os originais, que estes não 
contenham erros, falta de dados etc., pois 
que, na fase de revisão de provas, já não se 
admite inclusão de novo original e apenas 
se corrigirão erros de composição tipo- 
gráfica. 


10 As referências bibliográficas serão 


indicadas no texto por ordem numérica, 
usando algarismos árabes entre colchetes [], 
e a lista das mesmas, no final do artigo, 
obedecerá às seguintes normas: 


Artigos de revistas: - 
Nome(s) do(s) autor(es) e iniciais (ano de 
publicação entre parêntesis). Nome do 
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artigo, Nome da revista, número do 
volume, página(s). 


Livros: 

Nome(s) do(s) autor(es) e iniciais (ano de 
publicação entre parêntesis). Título, 
editora, página(s). 


Convem distinguir claramente entre 
referências e outra bibliografia que se 
recomende. 

Trabalhos em que a matéria apresentada 
não esteja referenciada não são, em 
rincípio, publicados. Este é um ponto 
undamental em que se insiste. 


11 Os quadros ou tabelas, numerados em 
romano, devem figurar no texto 


12 


dactilografado, podendo ser encimados por 
um título conciso. 


Os gráficos devem, de preferência, ser 
desenhados a tinta da china, em folhas 
separadas e de modo a estarem prontos 
para serem ERTEO com a maior 
nitidez e qualidade gráfica possível. 
Quaisquer letras ou algarismos que neles 
figurem deverão ser suficientemente 
grandes para permitir a sua redução que, 
em princípio, fica ao critério do editor. 

As folhas com os gráficos, figuras ou 
fotografias, deverão ser numeradas, 
indicando-se o trabalho a que dizem 
respeito. Além disso, no texto, indicar-se-á 
o local para a respectiva inclusão. 


NOTA: A TÉCNICA salienta a necessidade 


de 


os autores, quando for o caso, obterem em 


tempo as devidas autorizações de reprodução de 
excertos de outras obras, sendo esta matéria da 
inteira responsabilidade do autor. 


A TÉCNICA oferecerá cinco exemplares da 


Revista, por cada trabalho publicado, 
independentemente do número de autores. 
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F. E. Rebelo Simões 


Aspectos Estatísticos da Análise Espectral Analógica 


F. E. Rebelo Simões 


Departamento de Engenharia Electrotécnica e de Computadores do |.S.T. 
Centro de Análise e Processamento de Sinais da Universidade Técnica de Lisboa (I.N.I.€C.) 


RESUMO 

Apresenta-se um estudo estatístico do sinal de saída 
do analisador espectral analógico de banco de filtros, 
para sinais de entrada gaussianos estacionários. 
Define-se um estimador para o espectro dos sinais de 
entrada, cujas descentragem e variância são 
discutidas, assim como o erro quadrático médio. 

Calculam-se as expressões aproximadas dos 3º e 4º 
cumulantes do sinal de saída e as suas versões 
normalizadas, em função dos parâmetros "Largura de 
banda de ordem nº" e "tempo de integração de ordem 
nº. Faz-se em seguida a comparação daqueles 
cumulantes com os correspondentes das distribuições 
x2 e lognormal. São feitas ainda considerações de 
índole numérica relativas a vários cumulantes de 
ordem superior, que apontam para um afastamento da 
distribuição do sinal de saída e do consequente 
estimador espectral já referido, em relação à 
distribuição habitualmente proposta do y2, e maior 
proximidade da lognormal, situação esta resultante 
essencialmente da utilização dum integrador 
exponencial. 

Finalmente descreve-se brevemente o funciona-- 
mento do analisador heterodínico e mostra-se que, em 
determinadas condições para o filtro de banda, todo o 
estudo estatístico anteriormente feito é extensível a 
este analisador. 


Introdução 


Em inúmeros problemas de Engenharia, Geofísica, 
Fisiologia e Econometria surge o conceito de espectro de 
potência e daí a necessidade de determinar esta 
grandeza empiricamente; isto pode fazer-se em tempo 
real directamente a partir das medidas feitas no 
momento ou, em tempo diferido, a partir de registos 
efectuados previamente. 

Nesta última situação, se dispusermos por exemplo de 
um registo magnético, tanto podemos tratar os dados 
digitalmente - o que abre maiores possibilidades de 
processamento - como passar o registo por um analisador 
analógico, o que o coloca na primeira situação mais 
restritiva. 


| Na realidade, como referiremos adiante, estes aparelhos dispôem 
em geral de 2 ou 3 filtros fixos em lugar de 1, cada um com uma 
diferente largura de banda. 


CEC 


Original recebido para publicação em 6 de Dezembro de 1989 


ABSTRACT 


A statistical study is presented of the output of an 
analog spectral analyzer of the bank of filters type, 
with a stationary Gaussian input signal. The bias, 
variance and mean square error of the spectral 
estimator thus obtained are discussed. Some 
approximate expressions are derived for the 3rd and 
4th cumulants of the output signal and their 
normalized versions, as functions of the “nth order 
filter bandwidth" and the “nth order integration time" 
of the low pass filter. Then a comparison is made 
between these cumulants and the cumulants of the y2 
and lognormal distributions. Some numerical 
considerations are made regarding the nth order 
cumulants, which point to the output distribution 
being more close to the lognormal distribution than 
the commonly suggested y2 distribution. This situation 
stems primarily from the use of an exponential 
Integrator. 

Finally the heterodyne analyzer is briefly described 
and it is shown that, under certain specified conditions 
for the passband filters, all the statistical study 
previously made is still applicable to this analyzer. 


É desta situação mais restritiva, mas mais expedita e 
económica, que nos ocuparemos. 

Vamos considerar dois tipos de analisadores 
espectrais: o analisador de banco de filtros, em que a 
frequência central de cada filtro percorre o espectro a 
analisar - continuamente ou por saltos - e o analisador 
heterodínico, de filtro fixo, em que inversamente é o 
espectro que se desloca até a frequência a analisar 
coincidir com a frequência central fixa do filtro do 
aparelho 

Quanto aos sinais a analisar, sabemos que os de maior 
interesse em Análise Espectral Analógica são os sinais 
estacionários (do ponto de vista da potência ou valor 
quadrático médio). Estão nestas condições os sinais 
periódicos (determinísticos ou não), de espectro discreto, 
e os sinais aleatórios estacionários em sentido lato, de 
espectro continuo. 

Vamos restringir-nos a estes últimos, dada a sua 
maior importância e a maior delicadeza da sua análise. 
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G (ua) 


figura 1.1 Analisador espectral de banco de filtros. 


1. O Analisador Espectral de Banco de 
Filtros 


1.1 Princípio de funcionamento 


O esquema simplificado dum analisador espectral 
analógico não heterodínico é essencialmente o indicado 
na figura 1.1. 

O processamento do sinal implica portanto os três 
elementos principais seguintes: 

- um filtro de banda estreita, centrado na frequência F 

em que se pretende estudar o espectro, e que produz o 
sinal v(t) a partir de x(t): 


+ mo 


x(1) A(t— v) di (1.1) 


vulto) = x(t) *h(t)= | 
— 
em que h(t) é a resposta impulsional do filtro. 
A variação de F consegue-se à custa de um jogo de 
filtros ou dum filtro único de frequência central 
variável. 
- um sistema quadrático que realiza a função ou 
operação 


AO =vA0) (1.2) 


- um filtro passa-baixo, de resposta impulsional g(t), 
que fornece o valor 


+ 


y(u glt— x) dv (1.3) 


A= y(t) * g(t)= | 


E 
No caso normalizado tem-se 


+ 

| gludi=1 (1.4) 
E é] 

e, admitindo-se que gít) é causal, o filtro passa-baixo 

comporta-se como um integrador (ponderado) , que 

realiza a média ponderada do sinal y(t)=v(t) , ou seja, 

da potência instantânea do sinal filtrado v(t). 

A saída z(t) parece assim, intuitivamente, dar uma 
indicação da potência média por banda de frequência, ou 
seja da densidade espectral de potência, mais 
vulgarmente chamada espectro de potência, do sinal de 
entrada x(t). 

Notemos no entanto que esta perspectiva é temporal, 
visto ser baseada no conceito de potência média 
temporal. Ora os sinais (de entrada) em apreciação são 
aleatórios, constituindo um processo estocástico 
(estacionário), ao qual corresponde um espectro definido 
a partir do conceito de potência média estatística. 


Concomitantemente os sinais de saída z(t) são 
igualmente aleatórios, constituindo eles próprios um 
processo estocástico igualmente estacionário!?), 

No modelo que vamos adoptar existe equivalência 
entre as médias temporais usadas e as correspondentes 
médias estatísticas (ergodicidade) , tanto para os sinais 
de entrada como para os de saída. Como é de esperar, a 
ergodicidade da saída resulta da suposta existente para 
a entrada. 

Esta ergodicidade dos sinais de saída permite-nos 
comparar com fundamento as previsões teóricas 
(estatísticas) com os resultados obtidos ao longo do 
tempo, nomeadamente as flutuações do sinal. 


1.2 Estimador espectral Sy(0) 
1.2.1 Determinação de E( z(t)) 


Vamos precisar as afirmações anteriormente feitas 
quanto ao significado do sinal de saída z(t), estudando as 
suas propriedades estatísticas e relacionando-o com o 
espectro que se pretende analisar, de forma a deduzir daí 
um estimador espectral. 

Este estimador ficaria completamente caracterizado 
pela sua distribuição estatística. Dado porém que a 
determinação desta representa um problema bastante 
dificil, limitar-nos-emos a estudar a descentragem, a 
variância e o erro quadrático médio, dando ainda 
algumas indicações sobre distribuições que se 
aproximam da real. 

Consideremos então um sinal real x(t) pertencente a 
um processo estocástico estacionário (em sentido lato) , 
de densidade espectral S,(w), e vamos procurar 
relacionar z(t), através da sua média Ejz(t)), com o 
espectro S;(w) em questão. 

Se atendermos às expressões (1,1)(1.4), tem-se 


+m 
E (u“(1) kelt— vw) di= 


Eldd))= | 
is (1.5) 


+ 


=R (0) | gt-udi=R (0) 


— Éh 


em que R(:) é a (autojcorrelação do sinal filtrado v(t), 
dada por 


2 No caso geral para que zít) seja estacionário em sentido lato é 
necessário e suficiente que x(t) seja estacionário de 2º e da ordens, 
devido ao elemento quadrático. No caso particular de x(t) ser 
gaussiano a estacionariedade deste processo em sentido lato 
implica a sua estacionariedade de qualquer ordem, o mesmo se 
dizendo para o processo z(t) da saida. 
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R ()= Elutt) v(t+v) (1.6) 


e RO) é a potência daquele sinal. 

Ora, se designarmos por S(w) o espectro de v(t), o 
qual é, como sabemos, a transformada de Fourier (TF) de 
Rot), tem-se 


+ di + 

E(dD)= | S(w) — = | S (0) |H (co) po (1.7) 
— DD ” 2m — 2n 

onde H(qw) é a resposta de frequência do filtro de banda, 


com (cf. figura 1.4) 


+ 
H(w)= | he IS de (1.8) 


— DE 


figura 1.2 Resposta impulsional do filtro passa-baixo. 
g(t-6) 


O 


figura 1.3 ara do filtro passa-baixo sobre a potência instantânea à 
saída do filtro de banda, y(t) = vt). . 

A expressão (1.7) é absolutamente ei e não se pode 
simplificar, a menos que se conheçam exactamente os 
dois factores da integranda. 

Vamos admitir no que se segue que o filtro tem 
resolução suficiente para o espectro Sw) a analisar, isto 
é, que a banda é suficientemente estreita de forma a que 
Sdw) varie lentamente dentro da gama em que |H(w)| 
toma os valores mais elevados. 

Podemos precisar um pouco esta hipótese, definindo a 
banda equivalente de ruído do filtro 


F. E. Rebelo Simões 


[Hop 
: 2n (1.9) 
5, —em 


* (Hx) 
FTUELI 

Vemos que B, representa a largura de banda de um 
filtro rectangular de ganho |H(w)|maz que produzisse a 
mesma potência que o filtro real produz para uma 
entrada de ruído branco. (Fig. 1.5) 

Com esta definição a hipótese anterior pode 
escrever-se: o espectro S;(w) tem pequena variação à 
escala de Ba. 

Nestas condições (1.7) toma a forma aproximada 
. 9 (1.10) 

Elanj=2B,S (2nF) IH (co) 
x FUI 
em que F é a frequência central do filtro. 

A relação anterior sugere que utilizemos para 
estimador espectral a grandeza 


z(t) 
2 
2B 1H) 


Vemos assim que a estatística de Sw) é, aparte um 
factor constante, a mesma que a de z (t). 


S (w)= (1.11) 


IJ 


IH (oo)? 


47, 
Õ 27 F 


figura 1.4 Espectro de potência do sinal de entrada e função de 
transferência quadrática do filtro de banda. 


IH (wo)? 


Bro] 


Jo F 


figura 1.5 Largura de banda equivalente de ruido dum filtro de 
banda. A área do rectângulo é igual à área do gráfico de 
[Hi co 2). 
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1.2.2 Descentragem 

A relação (1.10) é, como vimos, uma relação 
aproximada. A escolha de (1.11) para estimador 
mostra-nos que existe um erro de descentragem ou 
enviesamento (“bias”), definido por 

a A , 
DS ()=ELS (w)—- S (0) 
resultante da variação de Sw) dentro da banda do 
filtro. 

Se esta variação for pequena em valor relativo, pode 
obter-se com boa aproximação uma expressão para a 
descentragem. Para esse efeito usa-se o desenvolvimento 
aproximado (3 termos) em torno da frequência F, 


(1.12) 


2 


— 


(1.13) 


S (W)=8 (2nF) + 4—F)S (2nF) + S (2nF) 
(em que as derivadas são em ordem a [=w/2n) 
e substitui-se na expressão (1.7) de Efz(t)). Obtém-se 


Eldn)=2S (2nF) | [H (co) [2 df + 
O 


LE 


+28 uh) | (-F Hay + 
) 


+28 (2nF) | |H(w)|2 df + 
Rio (1.14) 
+ S (2nF) | (-F | H (o) PP df + 
ni 0 


ui 


+28 (2nF)(F -n| (=P) | Hw)|É df + 
k Ú 0 0 


+S (QnF)(F py | |[H(w)|* df 
x Ú 0 


em que F,y é a frequência do eixo de simetria 
(aproximada) de |H(w)|2. Notemos que nos bancos de 
filtros reais de banda constante existe em geral aquela 
simetria, tendo-se neste caso Fyg=F. Nos bancos de 
filtros de banda de percentagem constante existe 
assimetria em relação a F, mas há uma simetria 
aproximada em relação a uma outra frequência F,. Para 
filtros de banda igual ou inferior a 1/3 de oitava (ex: 10%, 
3% ), Fo-F <F (por exemplo Fg-F =0,7%F no caso de 1/3 
de oitava, e Fy-F =0,01%F no caso de 3%). 

A suposta simetria de |H(w)|2 em relação a Fy permite 
anular os 2º e 5º termos de (1.14). Podemos assim 
escrever, após substituição na expressão de DSydw), 
para wo=2nF 


A | (a 
DS (2nh)=(F,—-P)S (2nh) + S (2nF)+ 


I (1.15) 


A S (QnF) 


| 2 
2B k IH (00) az 
em que [,, dado por 


= | ([-F 2] H (oo) Pdf (1.16) 
0 

representa o momento de inércia de |H(w)|2 em relação 

ao eixo f= Fo. 

Notemos que, dada a proximidade entre F e Fy, as 
duas primeiras parcelas de (1.15) representam com 
bastante exactidão a diferença S,;(2nF 6) - S2nF), que se 
pode em princípio desprezar nas condições admitidas 
para Sj(w) (No caso dos filtros de banda constante, 
simétricos, as duas parcelas são exactamente nulas). 

Fica assim finalmente 

É 
2B, | Ho) 
PmuIX 


Pode verificar-se que 


(1.17) 


DS (Qnh)= S (2nF) 
TI I 


2 3 
VE (1.18) 


0 K 

em que K é uma constante que depende da forma do 

filtro. Para o filtro rectangular, K=12 . Para filtros 

reais, K é um pouco menor (para um Butterworth de 

ordem 8, K= 10.2, para um “raised cosine”, K= 7.65). 
Vemos assim que 


I 


A 
DS (enF)=— B2S (2nF) (1.19) 
x 9K ho x 
A descentragem cresce com o quadrado da banda e é 
tanto menor quanto mais rectilineo for o gráfico do 
espectro dentro da banda do filtro. 


1.2.3 Variância de Sy (wo) 

Até aqui não se fez nenhuma restrição aos sinais x(t), 
admitindo-se apenas que são estacionários (sentido lato). 
No que se segue vamos supor que o processo x(t) é 
gaussiano estacionário, de média nula “, 

Da definição do estimador tem-se 


variz()) (1.20) 


PONS (20h )| = eme 
i 4B; |H (|) 
a muIx 

E necessário portanto determinar var(z(t)). 

Já vimos que Efz(t)) não depende do tempo, o que é 
uma das condições de estacionariedade lata de z(t). Na 
realidade, como se pode provar, dada a existência entre 
x(t) e z(t) de duas relações de linearidade e de uma 
relação quadrática, a estacionariedade lata de z(t) exige 
a estacionariedade lata e ainda a estacionariedade de 4º 
ordem de x(t). 

Nestas condições pode definir-se a covariância de z(t) 

C(y= Eiláy-— Elo] lát+ w)— Elxt+ U)])= 
* (1.21) 
=R (Eta) 


em que R/1) é a correlação de z(t). 


d Esta hipotese da média nula não tira nenhuma generalidade ao 
problema, dado que se existir uma média Efx(t/j=0,o filtro Híw|) 
corta a componente contínua, fornecendo um valor Elu(t)j)= 0. 
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Se calcularmos a TF de C/:) - admitindo que existe - 
obteremos a parte própria do espectro de potência de z(t), 
isto é, sem a função delta na origem correspondente à 
componente contínua Efz(t)). 

Vamos designar este espectro próprio por S, (w), uma 
vez que se identifica com o espectro de zp(t)=z(t)-Efz(t)). 
Tem-se então 


+ 00 
S (w)= | C (ve "a (1.22) 
Z0 z 


-+ mo 


S, (03) e“! do/2n (1.23) 


C(uy= | À 


Por outro lado o filtro passa-baixo, que constitui o 
integrador, introduz a relação 


S ()=|G(WJS (0) 
*0 dd 


— BE 


(1.24) 


em que S, (w) é por sua vez a parte própria do espectro de 
YU“ eG(w)éa TF de g(t). 
Como variz(t))=C0), vem finalmente 


-+- DO 


d | 
varizt)= | GÊ S (o) — (1.25) 
— DO + 2n 
Resta assim determinar S, (w), TF de Cjt). 
Ora 
C (O=EwÃO t+ VE AO) = 
yu =Etu (8 v U)— [Etu — (1.26) 


= Elv*(t) vAt+ v-Ri0) 


A primeira parcela de (1.26) é um momento de 4º 
ordem do processo v(t) que, resultando do processo x(t) 
através de uma filtragem linear é, como este, gaussiano. 
Por outro lado, como já referimos, Efv(t)/)=0. 

Ora, para variáveis multinormais v; de valor médio 
nulo pode obter-se a relação [1] 

Ejo uv,v = q 
Elo u,iEtu,v + Elo v,u,u + Elo vu Htuçu,) 


Se fizermos 


27) 


DU = v(t) 


UV; = V, = v(lt+ 1) 


resulta facilmente, de (1.26) e (1.27) 
Ci) =2RG) 
y vu 
que é uma relação fundamental. 
Podemos agora calcular S, (w) a partir do 2º teorema 
da convolução 


S (wuj=2S8 (w)"S (w)/2n= 
7 [T v 


(1.28) 


+ co 
=2 | Spth) S two —A) dh/ 2x 


e O 


(1.29) 


+ o 
=2[ S,tA) S tA+0) dh! 2x 


— CD 


por ser S;(w) uma função par. 


4 A estacionariedade de z(t) implica a de y(t) da mesma ordem. 


F. E. Rebelo Simões 


Para um filtro rectangular e Sw) uniforme dentro da 
banda do filtro, Sy 09) varia linearmente. Para o mesmo 
filtro e Sy(w) variando lentamente dentro da referida 
banda, verifica-se uma ligeira modificação em relação à 
situação anterior, apresentando S, (w) junto à origem 
(f=0) uma concavidade voltada para cima ou para baixo, 
conforme a concavidade de Sw) esteja maiorita- 
riamente voltada para cima ou para baixo, 
respectivamente, dentro da banda do filtro. 

Para Sw) constante dentro da banda e para filtros 
reais verifica-se que a situação anterior se acentua, com 
a concavidade voltada para baixo junto à origem. 
Finalmente o caso mais geral em que Sy(w) varia 
lentamente conduz a um reforço ou a uma atenuação 
desta situação (Fig. 1.6). 


6] F 


figura 1.6 Espectros de potência dos sinais v(t) e yo(U = y(t-E(y(tU). 
É esta grandeza Syd) que introduzida em (1.25), e 
tendo em atenção (1.20) permitirá calcular var(S (2n:F)) 
var! S (2nh))= 


1 E 2 e À 
——— | Gs (wdf (1.30) 
4 Jo 


48) |H(co)] = 
FULL 

Notemos que tanto var/z(t)J) como var(S(w)) 
dependem de F através do espectro S, (w). 

A variância de S42nF) pode normalizar-se pela média 
E/S d2nF)):; 


” van S (2nk) (1.31) 


SE 
[E(S (Qro)? 


Quanto ao erro quadrático médio tem-se, por 
definição 


2=E [8 (Eh) -S en) Ea 
ou, expresso na variância e na descentragem, 
“=van 8 (QnF)+IDS (Quo)? (1.33) 
Sob forma normalizada, tem-se 
2 DS (2nF) | | 
gt p | add | (1.34) 


-— omnes E e] e 
[EIS (2n)p? EIS (rh) 
Todas estas expressões são gerais, não se podendo 

explicitar melhor, a menos que se concretizem as 

grandezas que caracterizam o problema, e que são: 
- O espectro S,(w) a analisar 
- À função de transferência quadrática |H(w)|* dos 
filtros de banda 
- À função de transferência quadrática |G(w)|* do 
integrador. 
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1.2.4 Cálculo exacto de vartSx(w)) em 
condições ideais 


A título ilustrativo, e para fornecer ordens de 
grandeza, vamos considerar uma situação em que é 
possível o cálculo exacto da varfz(t)/. 

Trata-se do caso em que o filtro é rectangular (módulo 
de H), o integrador é exponencial no tempo (situação 
comum, com um circuito RC) e o espectro S,(q) se supõe 
constante dentro da banda do filtro. 

À situação é ilustrada pelas figuras 1.7 e 1.8 


O F 


figura 1.7 Função de transferência quadrática dum filtro rectangular. 


Syo (qo) 


<---4BSy(27F)=A 


So —A/2 


2F-B 2F 2F4B 


O B 


figura 1.8 Espectro de potência do sinal yy(t)=y(t)-Ely(t)), para o 
caso do filtro rectangular da figura 1.7. 


tendo-se 
|G (w)|*= ana (1.35) 
1 + (wRC) 
S, (w)= 48? (2nF)[B-—| (| |f|< B 
0 
as à 4 
28 (Qu) [B—|/—-25]) If-2F|<B (1.36) 
282 (2nF)[B-|[+2F|]  |f+2F|<B 
O fora destes intervalos 


Nestas condições a expressão (1.30) é facilmente 
integrável, conduzindo ao resultado 


"= 9 | [1:37 
vans. (Qnh)j=5" (275) [U(Nn) = V (Na) 


com 


l 
DO = UA 2 D+ PD O (1.38) 


| | 
V(O=I2g (0-2g (0+g  ,O+g, (O) €1:39) 
E | m+1 m-—l 
f(B=k arctg hé VA 
| ns | 
s(8=> log 1 +R2 E) (1.41) 
N=2BRC (109) 
Me: (1.43) 
RE 
B 


em que Fy é a frequência central do filtro, quase 
coincidente com F. 

Para diferentes filtros de banda de percentagem 
constante, desde 1% até 1/3 de oitava, obtêm-se, até à 3º 
casa NET Us TETE] constantes do a 


ra [oo o 


quadro 1.1 


Pode verificar-se facilmente que se tem 


lim NU(Nn)=1 (1.44) 
N— 
ni NV(Nn)=0 (1.45) 
é 3 2 | 
lim NºV(Nn)=- (1.46) 
N> [ 
U(0O)=4 
v(0)=2 (1.47) 
A ' 
lim var S (2nh); =2S (2nh) 


N— 


De (1.44/45) conclui-se que var($ 42nF)) tem como 

expressão assintótica (N— 00) 
A É a | 
var(8 (2nh)) = —Sº(2nh) a 
EE US & N T 

Vemos assim que, tomando este valor assintótico 
como o verdadeiro, se comete um erro (por excesso), que 
para N> 10 é inferior a 9%, para N >25 é inferior a 5% e 
para N> 100 é inferior a 1,4%. 

Nestas condições pode-se escrever 


A |] ; 
vartS (2nh)=(1 ii Sº(2nF) (1.49) 
: BT + 


em que 69 é o erro de aproximação assintótica no cálculo 
de variS (2nF)), com0<6,<1,esefez T=2RC. 


10 Técnica - 91/17 


A grandeza T - tempo equivalente de integração - é o 
inverso da largura de banda do integrador, como se pode 
verificar facilmente [c.f. 1.56]. 


(1.50) 


A variância normalizada 02 é dada por (1.31), ou seja 


2 / (1.51) 
P=(1 0) — | 
o“ =( 2 BT 

Sd ss (1.52) 
ass BT 


Note-se que na hipótese assumida de S(w) não variar 
dentro da banda, não existe descentragem, pelo que 
E2=o2, 


1.2.5 O caso geral dos filtros reais 


Neste caso não se pode fazer uma análise exacta. Em 
vez disso podemos determinar expressões aproximadas 
para var(S(2nF)), pondo apenas duas hipóteses 
simplificativas: 


A. O espectro S;/(w) varia lentamente à escala da 
banda do filtro 

B. A banda do filtro é muito superior à banda do 
integrador. 


A hipótese B corresponde a admitir que G(w) está 
muito concentrada junto à origem, em comparação com 
Syilw), o que equivale a supor que à escala das 
flutuações de y(t), o filtro passa-baixo tem um tempo de 
integração 
flutuações. 

Esta hipótese está ilustrada na figura 1.9. 


“ grande, permitindo reduzir aquelas 


figura 1.9 Espectro de potência do sinal yo(t) = y(t)-Ely(t)) para um 
filtro real e função de transferência quadratica do filtro 
passa-baixo. 
Podemos nestas condições tomar, no integral (1.25) 
Sy ot) como constante e igual a Sy,(0), sem cometer erro 
apreciável. 


Vem assim 
" no i 
vartz (p=, (0) | | G(o) | Pe dim 
, Es (1.53) 
+ + ; 
-2 | Sim) df | | G(w)|2 df 
Usando agora a hipótese A, vem 
(1.54) 


S (0) = 8 (2nh)| Ho) e 


donde, recorrendo à relação de Parseval para o filtro 
passa-baixo e à paridade de H(«), 


F. E. Rebelo Simões 


+ Do 
var (auj=482 (2h | g“(t) dt | | H(co) “df (1.55) 
0 


— o 


Se definirmos largura de banda B, do integrador como 


MP (1.56) 
= [G(co)| PES (cf 1.9) 
É ca G(0) 
e utilizarmos (1.10), vem finalmente 
vartto, Bs (1.57) 


[Etxop" B, 
em que B, - banda estatística do filtro - é definida como 


2 


| | Ho) |2 df 
e (1.58) 


| | H(w)|* df 
É] 
É evidente que da relação S(2nF)= Cxz(t) resulta 
vartS (2nh) var (0) . é” = 
EIS Qnrp teor B, BY, 
em que T. é o tempo equivalente de integração, inverso 
de Ba. 


(1.59) 


De notar que, no caso do filtro rectangular, 


B=B,=B: 
No caso dos filtros reais pode provar-se que se tem 
B,> Bh. 


Por exemplo nos filtros dos analisadores Brúel e Kjaer 
tem-se 


B > B 
ro 3 h 
A expressão (1.59) mostra que a hipótese B pode ser 
escrita equivalentemente 
B «B ou ainda 
Ea (1.60) 
0*<] 
Na realidade a expressão (1.59) é, como vimos, uma 
expressão assintótica, devendo escrever-se com mais 
rigor 


var S, (2nF)) B 


=(1-0)-É=(1-0)— 
CEB < r. 


assis. aa (1.61) 
[EIS (Qi? 
Aliás a análise que acabámos de fazer contempla o 
caso do filtro rectangular estudado anteriormente, 
bastando pôr em (1.61) B,=B. O estudo exacto feito 
antes permite-nos, na realidade, determinar o valor de 
Os em função de BT. 
Voltemos à expressão do erro quadrático médio 
normalizado. 


à No caso do integrador exponencial T=2RC, como ja vimos. No 
caso geral Té por definição e inverso da largura de banda B,. 


Técnica - 91/1 11 


Aspectos Estatísticos da Análise Espectral Analógica 


Da expressão (1.34) de €2 conclui-se que para tornar 
esta grandeza arbitrariamente pequena é necessário 
actuar nas duas parcelas - a da variância e a da 
descentragem. 

Para reduzir a primeira é preciso aumentar 
indefinidamente o produto B,T.. 

Para reduzir a segunda é preciso inversamente 
diminuir B, (cf 1.19), e portanto B,. 

É óbvio que estas condições implicam 

B >0 

E.) 
a —> 00 (1.62) 
BT >» 

E E 
o que se consegue fazendo decrescer mais depressa 1/T, 
ou seja B,, do que B,. Podemos portanto concluir que, se 
fizermos variar B, e T, de acordo com as condições 


anteriores, o erro (relativo) do estimador $;(2nF) pode 
tornar-se tão pequeno quanto quisermos. 


Isto é 
, ta 
lim EIS (2nF))=0 
B T P —+ E (1.63) 
8 E e 
B > 
B 
ou ainda 
(1.64) 


Lim. 8 (en) O=S (2nF) 
Isto significa que o estimador S,(2nF) é consistente. 


1.83. 
1.3.1 


A determinação exacta da estatística de z(t) - ou de 
Sw) - é um problema naturalmente mais difícil do que 
a da variância, que já considerámos, sendo dum modo 
geral inviável analiticamente. 

A literatura da especialidade apresenta duas atitudes 
diferentes nesta matéria: 

- uma simplista , mais corrente, em que se aproxima a 
distribuição por uma distribuição em yí - mais 
precisamente - Ar,? - em que o factor À e o número v 
de graus de liberdade são determinados de forma a 
fornecer a mesma média e a mesma variância que 
z(t). 

- outra, muito pouco frequente, em que se determina 
numericamente - dadas certas hipóteses relativas ao 
filtro ou ao integrador - a distribuição, através por 
exemplo do cálculo numérico da transformada de 
Laplace da função geradora de momentos (ver por 
ex: [2)). 

Vamos tomar aqui uma posição um pouco diferente - 
vamos, embora sem o rigor numérico da segunda 
atitude, procurar ir mais além do que se faz na primeira 
atitude, isto é, vamos determinar aproximadamente 
alguns momentos de ordem superior de z(t) e 
compará-los com os correspondentes da distribuição y”, a 
fim de concluir da validade de usar esta distribuição 
como modelo. 


Outros parâmetros da distribuição 
de z(t) 


6 Lim. - limite em média quadrática (limit in the mean). 


1.3.2 - 0 3º momento 
L3.21. 


Na realidade vamos determinar o 3º momento central 
de z(t), que coincide com o 3º cumulante (precisaremos 
adiante este conceito). 

Seguiremos um caminho paralelo ao seguido na 
determinação da var/z/, que é um 2º momento central. 

Assim vamos obter este momento kg” a partir do 
espectro (próprio) de 3º ordem de z(t) - Transformada de 
Fourier dupla da função cumulante Ky“(ry,19) - dado que 
ky = Ks(0,0). O espectro de 3º ordem de z(t) relaciona-se 
com o de y(t) através de uma função de transferência de 
3º ordem. Será pois necessário determinar o espectro 
Sa (w1,w9), o que se faz por TF dupla do 3º momento 
central (cumulante) Ka'(ry,19). 

Como y=v2, o referido momento é na verdade um 
momento de 6º ordem do processo gaussiano de média 
nula v(t). Graças às propriedades dos processos 
gaussianos aquela grandeza calcula-se fácilmente, 
assim se completando a operação. Veremos a seguir em. 
pormenor como se procede. 

Vejamos então. Para a 3º ordem a função cumulante 
identifica-se com a covariância de 3º ordem, tendo-se 
para z(t) estacionário ” 

Kil)= Ci t)= 
lixo — ElxoiAt+ 1) — Elo lat+ vo) — El) 


Pode definir-se o espectro de 3º ordem ou bispectro de 
z(t) 


(1.65) 


+ 0 ii 
is aj tado 1 
Sg (0,05) =)) esp * 


que com K3'(L4,19) constitui um par de Fourier 


Lj + iRgto) E, dr, (166) 
E 


Két 


gp! 


Ir +05 (to ,L, FO) 
==) “efitajage E às do, duo, (1.67) 
4 
Por outro lado, e em relação às grandezas 
correspondentes do sinal y(t), tem-se, de (1.3) 
Si (000) =S) (0,05) T (00,009) (1.68) 
em que Td(wy,wo) é a função de transferência de 3º 
ordem do passa-baixo, que é dada por 


(1.69) 


A grandeza pretendida k'; que estende para a 3º 
ordem o conceito de variância, obtém-se a partir da 
função cumulante com 1 =19=0. 

Portanto 


Toya) = G(w,) G (005) Gl —wW, — Wo) 


E SEN. A Eae. 
Fi=k5001=— | | SYoyapT oo du do, (1.70) 


T Ver nota (2). 
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O problema que se põe então é a determinação de 
Sy(wr1,w2), ou seja, do bispectro (próprio)” do sinal y(t). 
Para esse efeito vamos calcular Kg(try,19). 


Se fizermos 
ytb= 3, 
poa di (1:71) 
y(+)=3, 


e atendermos a que a função cumulante de 3º ordem se 
identifica com a função momento central, vem 
fácilmente 


K 13209 (1.72) 


Eb) (Eb Eb) +2Eb, Eyes) 


em que ( )y representa todas as permutações possíveis 
sem repetição. 

O primeiro termo do segundo membro desta expressão 
é o momento de 6º ordem Efv2(t) v2(t +11) v2(t+19)). Para 
o cálculo deste momento utilizamos a relação geral, 
válida para o caso gaussiano [1] de média nula 


Elo v,vav VV g= (Elo uphEtv o Elo uh : (1.73) 

Quanto ao segundo termo, é constituído por 3 parcelas 
do tipo Ef v2(t) | Ef v2(t+ry) v2(t+ 19) | produto dum 
momento de 2º ordem por um de 4º ordem [cf (1.27) ]. 

Efectuando os cálculos e reduzindo os termos 
semelhantes, chega-se sem dificuldade ao resultado 

K; (1 ty)= 8 Rot) Rio) R (i— L) (1.74) 

que é uma primeira generalização de (1.28). 
1.3.2.2 Bispectro (próprio) de yí(t). 
Cálculo de k73 

Trata-se como dissemos da TF dupla da expressão 
anterior. 

Tem-se assim 
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integrador. Admitimos assim que concentração de 
Talwr,w2) junto à origem é muito superior à de 
Sgy (w1,w9), e vamos tomar em (1.70) S3(0,0) em lugar de 
S3'(w1,09). É de esperar que o erro cometido, por excesso, 
para uma mesma relação B,/Ba, seja superior ao da 2º 
ordem, dada a subida de dimensão (basta comparar, para 

o caso do filtro rectangular, as concentrações de T, e Sy' 
ao longo dos eixos principais e da bissectriz wj = w9). 

Virá então 
poe (1.77) 


E | e 
k4=( —8,) Sa(0,0) | [o E kom) dw dw,/ 4n 


+ 00 + 0 
Ri 3 | | 2 
=(1-0,)8 | = af | IN T om) do dw/4n 


em que 05 é o erro cometido (0 <B4< 1) (9), 
Ora 


+ ao 
| Si(w) df= 28º Enf | Hw)º df (1.78) 
a U 


Por outro lado 


Ra 
eme | n TÃO 0) do dw, 


+ mo 
- | g“(t) dt 


como se pode verificar, exprimindo G( --wg) em g(t). 
Finalmente 


kg=(1-0,)16 S(2nF) | [H(o))É df | g3 (8) de (1:80) 
0 0 


para o caso real em que o integrador é causal. 
À expressão (1.80) com 03 =0 representa a expressão 


Sa (0 pus) = assintótica de k'3, quando B,/B, tende para zero. 
+ =Jtw LL +Hwoto) (1:75) 1,3.3 O 4º momento 
=8| [Ra pRpR pe dito 1.3.3.1. 


Exprimindo por exemplo Ry(tg-1;) na sua TF, Sw), 
resulta finalmente 


= | a (1.76) 
Sal iso) =8| PA 5 Ah+ 04) S A+w, li 
que é a primeira extensão da expressão (1.29). 

Para o cálculo de ky' teremos de recorrer a (1.70). 

Tal como na 2º ordem, esta expressão só se pode 
determinar com rigor se se conhecer exactamente S,y(), 
ou seja, Sdw) e H(w), podendo no entanto calcular-se 
aproximadamente se se fizer o mesmo tipo de 
simplificações que se adoptaram para a 2º ordem. 

Assim, vamos admitir que Sy(w) varia lentamente 
dentro da banda do filtro e que esta é muito superior à do 


8 Obispectro calculado a partir da funç ão de corre lação de 3º ordem 
do sinal y(t), de média diferente de zero, contém funções deltas. E 
nesse sentido um bispectro impróprio. 


Ao contrário do que sucede nas ordens inferiores, o 
momento central é distinto do cumulante. Sigamos um 
caminho análogo aos anteriores e comecemos por 
calcular a função cumulante de 4º ordem de y(t). 

Esta função é igual à covariância múltiplalll das 
variáveis 


y,= (8 
Yy= ytU+) 
(1.81) 
Y,=90+19) 
9 IU+T) 


8 Notar que, para uma mesma relação B,/B, se deve Ler 
o<B,<bB,<l 
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Por definição esta grandeza obtém-se por recorrência, 
por subtracção sucessiva de todas as covariâncias 
múltiplas de ordem inferior, a partir do momento de 4º 
ordem. 


Kg) = ED yoga) A O) Goa) | 
aa = - (1.82 
[E 90,99) Kad) o, [KO KO) Ko a04) e ( ) 


=K OK yo) K (9) KO) 
com 


” tu RÍ = (1.83) 
KO)=Ety) e KO I)=R(0,%o) 


e em que como sempre as chavetas indexadas indicam 
todas as permutações que conduzem a resultados 
distintos. 

Usando os resultados da 2º ordem e da 3º ordem 
obtém-se 


EO) = E a) Eb Et 

- Eb MEG JEb EO) 

O primeiro termo do 2º membro é o momento de 8º 
ordem Ef v2(t) v2(t+1;) v2(t+19) vê(t+1g) |. Para o 


cálculo destes momentos usamos a relação geral, válida 
para variáveis gaussianas de média nulal1] 


(1.84) 


Elo o,0,0 0 sorv,b=( Elo voo o JE(o o JE(o od] Re -B5) 
constituído pelas 3x 5x7= 105 permutações possiveis 
dos índices 1 a 8 sem repetir parcelas. 

Os restantes termos fazem intervir momentos de v(t) 
de 2º, 4º e 6º ordens, já estudados. Substituindo as 
expressões de todos estes momentos em (1.84) e 
reduzindo os termos semelhantes, obtém-se, após 
simplificação 


K'(r 


4 porta) = o IR tp R (o) E Mg— L|) Ro, a Lo) + 


+RO)R MR (DR (1, —1)+ (1.86) 
+RO)RA)R (Ri) 
1.3.3.2 Trispectro (próprio) de vít). 
Cálculo de k”4 


Trata-se agora da TF tripla da grandeza anterior. 
Exprimindo Rj(tg-19), Rutg-t1), Roltg-19) no espectro 
Sul), resulta finalmente 
y im | 
S 4 (009, 00,)= (1.87) 


16 lo ww) + (mo 0,09) +] (mo, 49,)] 


em que 
Ea (00,,0,,0,) — 
=| SAS A+ro)S Aro +w,). (1.88) 
ão v | 2 
dh 
SA+ro to, +w,) us 
14 


e que é uma extensão natural do integral de (1.76). 
Notar a relação de simetria 


a (1.89) 
Lilo 0,04) =D (000) 
Procedendo agora como para a 3º ordem, vem 
15 Em 
SO Do 4) = (1.90) 


qa tr , FIN 
= Sql, Wo, 094) Ao, W,, w,) 
com 


E AO, Wo (9) — 


(1.91) 
= U (0) G (0,5) G (0,5) G (— wo, —W,—,) 


função de transferência de 4º ordem do integrador. 
E finalmente 


k4 =K410,0,0 = (1.92) 


I [E a] y E | 
E 2 4 | | Soyogag To ag do do sda, 
(Zn) 


que se deve comparar com (1.70). 

Para o cálculo exacto desta grandeza põe-se as 
mesmas questões que pusemos nas ordens anteriores. 
Resta-nos assim igualmente o recurso ao cálculo 
aproximado, duma forma semelhante ao que já foi feito, 
e que - vê-se intuitivamente - vai conduzir a um erro por 
excesso, superior aos das ordens inferiores. Tem-se assim 


T tw, ww )do, dw, du, 


E Y E E 1 23 2 3 
É = 010,89 10,0,0) | | pot 


3 
Hu (1.93) 


Como 


| S(w)df=28º cen | Htw)Sdar (1:94) 
—w 'Q 
e por outro lado se tem 


| 
(211 é 


| | om (09, Wo, wo) dw, do, dw, = 
(1.95) 


= | sds 


resulta finalmente 


Pres ; +d " é ú 1.96 
há =(1-0,)965S (2nh) | [H (qo) | ar | g*(1) as | 
0 U 


para um integrador causal. 

Chegámos assim a uma expressão homóloga de (1.80), 
sendo de esperar que se tenha 0, >03. 

Note-se que 05, 03, U4 podem fazer-se arbitrariamente 
pequenas, à custa da pequenez do valor B9/Ba de que 
dependem, ou seja, da pequenez de o2 (ou de €2, se 
desprezarmos a descentragem), ficando as expressões de 
k'o, R'3€ k 4, com 0,=0, como expressões assintóticas, 
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1.3.4. Resumo dos cumulantes estudados. 
Generalização. 


| S2(9nb) [ (u q 1.55 
hi =4(1-0,) 82 (2nF) | ro) | edi 59 
0 0 


| 9 ç DM im (1.80) 
ko = 16(1-0)5 (2nh) | Hof | g“udt 


O 


li ii (1.96) 
i= 96-09 8!En 1 | Ho Par | g“ (di 
* 0 0 


com 


0 


As expressões anteriores são generalizáveis, embora 
com erro crescente com n, podendo obter-se [3] 


lim 
Bob 
É 


pé =9"(n— DIS Hmh) | Hc” df | grado ) (1.97) 
R x Ú () ft 
Notar que também aqui 
lim O =0 
9 A 
o —l 
e que se tem formalmente k*; = Efz(t)), representando 0; 


a descentragem normalizada (neste caso podendo ser 
0,<0). Fazendo 


R - | Ho)" df (1.98) 
FL 0 
2 =| g"(t) di iade: 
0 
(1.97) escreve-se 
k=2n-1!S2nF)P Q (1-0) (1.100) 
FL A no R n 


ou, sob forma normalizada, e desprezando o erro de 
descentragem 6,. 


k, Ei n 
E = =p =D.) (1.101) 
(k?) PQ, 


A expressão anterior pode tomar outro aspecto, se 
introduzirmos os conceitos de largura de banda e de 
tempo de integração de ordem n, respectivamente B, e 
Ta 

Assim, se definirmos estas grandezas , para n=2, 
através das relações 


e | =" (1.102) 
B, > Pp (10) 


n 


10 É imediato queB, =B,eT,=T, 
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A (1.103) 
nm q, 
resulta 
(int= 1)! 
é=(1-0)— (1.104) 
Ro jr 
rn mn 


que se deve comparar com a expressão homóloga de 2º 
ordem (1.61). 


1.4 
1.4.1. 


Conforme dissemos no início, vamos começar por 
seguir a via já referida em 1.3.1., de comparar os 
cumulantes calculados aproximadamente (n=2, 3, 4) 
com os cumulantes da distribuição xy? com a mesma 
média e a mesma variância, ou seja, com a mesma 
variância normalizada. 

Para 2z= Ar, tem-se de um modo geral [4] 


Comparação com a distribuição y2 


+ 


EX =9"A"y(n—1) (1,40) 
n 
donde, sob forma normalizada 
2 9 (1.106) 
a tn — 1) 
a 
Em particular 
é q 4/4 
O q (1.107) 


pelo que se pode escrever, explicitando esta grandeza, 


2 = 
c, =(n- 10? 
Comparando agora (1.104) e (1,108), e tendo em 
atenção que a igualdade de 0º nas duas distribuições 
implica 


(1.108) 


o BR l | 
0*=-=(1-0,) (1.109) 
v Bo! 
resulta 
Z E sm n-1] 
A quad (1.110) 
Ê Ea (= )0=4 Bl, 
! RE E 
ou 
2 
E 
ta FA [:17t) 
pt 
em que f,(09), dado por 
1-0 
f (0,)= n (1.112) 
n 2 n-l 
1-0) 


é o factor de correcção em relação à expressão assintótica 
de k',/k%, (próximo da unidade para n=3 en=4e para 
021) e r, é um coeficiente que se pode decompor em 
dois factores p, € qn, devidos respectivamente ao filtro e 
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ao integrador 


a) (1.113) 
B 


| il 
(=) (1.114) 
a=(7") 
n 


Como se vê, este factor r, é o único responsável, na 
situação assintótica fn(09)=1, pelo maior ou menor 
afastamento da distribuição de z(t) em relação ao yº. 

Quanto ao número de graus de liberdade v da 
distribuição xy? com a qual estamos a comparar a 
distribuição de z(t) tem-se, de (1.109) 


(1.115) 


Consideremos agora o caso particular do conjunto 
filtro rectangular (na frequência) + integrador 
rectangular (no tempo) : tem-se B,=B, T,=T, pelo que 
hzn/hX, =fn(09), isto é, a distribuição de z(t) é 
amtintnticamente x2, com v=2BT graus de liberdade. 
Repare-se porém que, à medida que o? aumenta (e com 
ela 69), a distribuição se vai afastando do y2 graças ao 
factor fn(09). 

Este factor é difícil de determinar no caso geral, para 
qualquer filtro e para qualquer integrador. Tanto o 
factor 1-0, como o factor (1-09)n-1 diminuem quando n 
aumenta, para um mesmo B5/By : o primeiro porque o 
erro de aproximação usada é cada vez mais grosseiro, o 
segundo porque 1-0, <1. Resulta deste modo uma certa 
compensação entre as duas diminuições. 

Com a finalidade de determinar a ordem de grandeza 
deste factor fn(03) foi estudado exactamente um exemplo 
com filtro rectangular e integrador não causal com 
resposta de frequência exponencial, obtendo-se para 
fnl(09) valores como os que a seguir se exemplificam no 
ETs 1.2. 


o2w|BT [1.02] f3 | fa | fo | fo | fro| fis] fzo 
Tolo eso 557 [38 [pes [em [on mo[3o 
E lofesf fm [miojo mo 
5 o [aos | one on | e oo | sa 0a [01 


quadro 1.2 

Para filtros reais f,(09) diminui em comparação com o 
caso rectangular. No trabalho já referido, dada a 
dificuldade de um cálculo exacto foram feitas certas 
conjecturas sobre a natureza do factor fn(09), obtendo-se 
nesta hipótese uma redução moderada daquele factor, 
para um filtro de 1/3 de oitava (analisador Brúel & Kjaer 
2113) e para um filtro de 300 Hz, menos rectangular que 
o anterior (Brúãel & Kjaer 2010). 

Verifica-se que para o? suficientemente pequeno 
(02 <0.02) se podem obter erros pequenos até ordens 
elevadas. No caso presente em que apenas comparamos 
cumulantes de 3º e 4º ordens, o erro em causa parece ser 
da ordem de 2% para 02 = 0.025, (0= 0.15) sendo portanto 
perfeitamente aceitável o uso da expressão assintótica. 


Foram feitas determinações de B, e pa para vários 
filtros, de bandas nominais entre 3 e 300 Hz e ainda de 
1/3 de oitava centrada em 1000 Hz, dos analisadores 
Brúel & Kjaer já referidos. Para n=3 e n=4 
obtiveram-se os resultados do quadro 1.3 


RESENDE 
DEDEDED 
[epa 
MEMES 


quadro 1.3 


Por outro lado, para o integrador RC usado naqueles 
analisadores tem-se 


(1.116) 


T=n* É RC 
donde 
qu=a (1.117) 
= 
resultando 
4 
= UF 


O valor global rn no caso dos filtros ex: ca ii é 


[ 55: [30H] sonavo 
RICO] 


[e ass [ren [me 


“quadro? 4 
Vemos que a influência do integrador RC é muito 
mais importante que a dos filtros reais no que respeita 
ao afastamento em relação ao y2. Este facto acentua-se 
drasticamente com a ordem n. 
Por exemplo para o filtro de 300 Hz, obtém-se 


Pio = 4.30 


Pis = 13.03 


is = 1092 


qro = 51.2 |q20 = 26214 


quadro 1.5 
Os cumulantes de 3º e 4º ordens permitem calcular 
dois parâmetros relevantes : o coeficiente de assimetria e 
a kurtosis. 
Para esse efeito definem-se os parâmetros 8; e Bo. 


2 

(ka) (1.118) 
a 3 
o) 

j = k, (1.119) 
2 


= = 4 * i 
O coeficiente de assimetria - "skewness" - é definido 
como o quociente entre a diferença média menos moda e 
o desvio padrão, e mede a assimetria relativa da 
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distribuição. Verifica-se que para uma larga família de 
distribuições, esta grandeza é dada pela expressão 


SE Gra 


Dada a relativa complexidade da expressão, é mais 
usual utilizar simplesmente a expressão 


Sk,=VB, 

que também dá uma boa indicação da assimetria, 

Quanto à kurtosis pode ser definida pelo próprio 
parâmetro 59 e mede o maior ou menor achatamento do 
pico da função densidade de probabilidade em relação à 
distribuição normal : mais achatada se 59<0, mais 
aguçada se 59>0. 

Calculemos então 5; ef» no caso y2. 

Ora de (1.105) tira-se 


[5] (1.120) 


(1.121) 


8 
R,=—=4e” 
W 
pelo que 
a (1.122) 
Sk=[B, =20 
g É (1.123) 
2 w 


Vejamos agora as expressões de V$ e Bopara a saída 
z(t) do analisador. Recorrendo a (1.104), dado que as 
formas adimensionais de Lj e 52 permitem usar 
cumulantes normalizados, e utilizando (1.112), vem 


ns | (1.124) 
/B, =21,09r,o 
(1.125) 
= 2 
B;= 6f Ar 4º 
Concluímos assim que 
(Sk)f | 
—a Ig, (1.126) 
(Sky 
(Bo) 
: =[ Or, (1.127) 


(By 


Para 02<0.1, fy(09) e f4(05) são muito próximos da 


unidade : dentro das conjecturas feitas são superiores a 


0.96 e 0.93 respectivamente, para os filtros habituais. 
Nestas condições os valores de (1.126) e (1.127) são 
apenas ligeiramente inferiores aos do quadro 1.4, 0 que 
significa que o coeficiente de assimetria e a kurtosis são 
cerca de 1.4 vezes e cerca de 2.3 a 2.4 vezes 
respectivamente os correspondentes à distribuição y2. 

Por outras palavras, tanto a assimetria da função 
densidade de probabilidade de z(t), como o seu 
aguçamento em torno da moda são mais acentuados do 
que na distribuição y2. 
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Concluímos assim, quer pelos parâmetros anteriores 
derivados dos cumulantes de 3º e 4º ordens, quer 
principalmente pelos cumulantes de ordem superior (ver 
Q 1.5), que a distribuição parece afastar-se da 
distribuição yº, mesmo se tivermos em conta as 
correcções que se devem introduzir através do factor 
fnl09) devido aos erros de aproximação cometidos no 
cálculo de kz,. 

1.5 Comparação com a distribuição 
Lognormal 


A função de distribuição lognormal tem um 
andamento semelhante ao do y2, apresentando no 
entanto momentos e cumulantes que crescem muito 
mais depressa com a ordem do que nesta distribuição. 

Dados os valores elevados já referidos encontrados 
para k,º nas ordens superiores, é natural comparar as 2 
distribuições (z e lognormal), pelo menos para alguns 
valores de n, a fim de ver se este modelo é viável. 

Embora alguns estudos anteriores |3|] mostrem a 
proximidade dos cumulantes das 2 distribuições para 
valores de n até ordens relativamente elevadas (n= 20), 
vamos limitar-nos a considerar, como fizemos com o yº, 0 
3ºeo 4º cumulantes. 

Assim, para a distribuição lognormal, a partir da 
expressão do momento de ordem n, pode obter-se 


k5" =aºo” 
RN =a'o%3 +02) roi) 
2a | 
RO =a*o*16+150"+60" +09) 
em quea=Ejz). 
Daqui resulta 
c; 2 (1.129) 
J E 
O sabor 
C 3 
3 
c, 6 (1.130) 
4 
"EN 2 Ra PU 
º l6+150"+60 +o 


“4 


Os valores do quadro 1.4, não considerando o factor 
[n(09) passarão a ser, parao2=0.1 e 02=0.05 
02=0.1 02 = 0.05 


AEEIEDEE 


quadro 1.6 


muito mais próximos da unidade do que no caso y2. 
Da mesma forma se obtém 


: (1.131) 
Sk =/B, =(3+02 
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h (t) 


FILTRO DE 
BANDA 


a cos (upt+a) 
figura 2.1 Analisador espectral heterodinico 


| 1.132 
B,=(16+150"+60!+09%o” 
Daqui se conclui 
(Sh)N J J 
(Bo o (1.134) 
(BN 48 dra 


Nas mesmas condições em que comparámos as 
distribuições de z e x2, ou seja, 02 <0.1, verifica-se que 
tanto a assimetria como o aguçamento em torno da moda 
da distribuição de z(t) são apenas ligeiramente inferiores 
aos da distribuição lognormal. 

Se acrescentarmos a isto o facto de que r, decresce 
com n muito mais lentamente do que r, crescia com n - 
por exemplo, e usando valores assintóticos de r, , tem-se 
rio =0.8, ris =0.6, rogo =0.5, para filtros de 300 Hz 
(comparar com os valores r10, r15, r99 do Q 1.5) - mesmo 
tendo em conta a redução imposta por f,(65), concluímos 
que a distribuição procurada de z(t) é muito mais 
próxima de uma lognormal do que de uma y2. 


2, O Analisador Heterodínico 


2.1 Introdução. Estacionariedade 

O esquema simplificado dum analisador espectral 
heterodínico está representado na fig 2.1 e corresponde à 
introdução de um modulador de amplitude no esquema 
da fig 1.1. 

As relações entre v(t), y(t) e z(t) mantêm-se [(1.2), 
(1.3)], sendo (1.1) substituída por 


Da 


figura 2.2 Espectro de potência do sinal modulado u(t) 


g (t) 


FILTRO 
PASSA BAIXO 


(2.1) 


ult) = ult) * h(t) 
e acrescentada a relação 


ult) = axtt) costa t + à) Ra:8] 

Para sinais de entrada estocásticos estacionários, O 
modelo adoptado para o modulador possui uma fase a 
aleatória com distribuição uniforme num intervalo de 
largura 2n, o que é de resto intuitivo. Este modelo 
permite proporcionar aos sinais a juzante, u(t), v(t), y(t) 
e z(t) um carácter igualmente estacionário. Com efeito, 
pode provar-se [6] que a estacionariedade lata de z(t) 
implica a estacionariedade de 2a e 4a ordens de x(t), tal 
como no caso não heterodínico já tratado. 


2.2  Estimador espectral S(c)) 
2.2.1 Determinação de E(z(t)) 


De (2.2) deduz-se imediatamente, e para x(t) 
estacionário em sentido lato, 


9 
R ()=—R (tcos wL (2.3) 
LÊ 2 X 
o que implica, por simples TF 
2 
(2.4) 


a 
S (w)=—IS (w-w)+S (w+ w,)| 
u 4 x O X O 


espectro este representado na fig 2.2 

No analisador heterodínico wp( =2nF ) é variável com 
a frequência F em que se pretende medir o espectro de 
x(t), tendo-se 
A (2.5) 
F = F E F 
em que F; é uma frequência fixa (a frequência central do 
filtro fixo de banda estreita). Se F; for muito superior à 
maior frequência em que S,(w) toma valores apreciáveis 
(o mesmo se dizendo portanto de F,) as duas parcelas de 
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(2.4) são disjuntas(1D), Por outro lado, se S,(w) varia 
lentamente à escala da banda do filtro, ter-se-à 


2 


S ()= TS (PH (2.6) 


em que F é a distância entre Fye F. 

Para o cálculo de Efz(t)) recorre-se de novo a (1.2) e 
(1.3), vindo para o integrador normalizado (1.4) uma 
expressão homóloga de (1.10) 


2 


as À a LED é NB (2.7) 

Eixo) = R (0) = : S (2nF)B O 

que se identificará com (1.10) se fizermos a =2. 

Daqui resulta para estimador de S;(w) 
A 2xt 
É (ini) = om (2.8) 
ú a“B |H(w)|? 
h “max 


Dada a identidade das expressões de Efz(t)) nos dois 
tipos de analisador é obvio que todas as considerações 
feitas sobre a descentragem se aplicam igualmente nos 
dois casos. 


2.2.2 Variância de Sw) 


Vimos no caso do analisador não-heterodínico que se 
tinha, para x(t) gaussiano, e portanto para vít) 
gaussiano, a relação 


C()=2Rjy) 

No caso heterodínico, quando x(t) é gaussiano, uít) 
não o é, devido à modulação, o mesmo sucedendo no caso 
geral a vít), após a filtragem linear, pelo que a relação 
anterior, base da análise da variância de z(t), deixa de 
ser válida. 

Veremos no entanto que, sob certas condições para o 
filtro - existência de faixa contínua de rejeição total(12) - 
é possível que a expressão anterior se mantenha, e com 
ela a validade da teoria anteriormente feita. 

Vamos então calcular Cr) a partir da sua definição, 
exprimindo y(t) em x(t). 


Cy = Ely(bylt+u) — Ely(OE(y(t+ v)) 


(1.28) = (2.9) 


(2.10) 
=Eju“(gvt+ v) — EjuAOE(u “t+ 1) 
em que, de (2.1) e (2.2) se tem 
+= DD 
db=a | hq) x(t— pu) cos lo lt- pl +aldu 82:11) 
Se fizermos 

t=t-p, tg=t+1—-n,=Ê, (2.12) 

to=t-Ho Ethel St 


1i Na prática, para espectros de audio-frequência, usa-se por 
exemplo F, = 1.2 MHz, com banda de 10 KHz. Modula-se a seguir 
em F= 1.23 MHz e filtra-se a 30 KHz, com bandas entre 300 Hz e 
1000 Hz. Deste modo evitam-se sobreposições de espectros. Pode 
prosseguir-se este processo para bandas mais estreitas (3 a 
100 Hz) modulando de novo com F,=30.75 KHz e centrando os 
filtros finos em Fl =750 Hz (ezemplo de Brúel & Kjaer 2010). 


12 Fisicamente irrealizável, mas aproximável na prática. 
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xt )=x 
otta=b :=12,3,4 (2.13) 
cos B. =c, 
vt )=x", 
ont ta=p i=1,2 (2.14) 
cos B, =€. 


resulta, tendo em atenção a independência suposta 
existente entre os processos x(t) e cos(wgt + a), 
+o d 


Ciy=at| [| nu) E| [15,Je| [Ie] - mad 
-s[Me efe o(s efe) ]to, 


Como já vimos [cf. (1.27)] 


4 
2.16 
E| | | x) = [Eta ee), (2.16) 
Por outro lado 
4 
E| [ej=c,+0,+c, aim 
2 3 [É 
| | | E] [| e =0,+06, 
com 
8C =cosw(21+pn +n,—u,—H))] 
8C,=cosiw (Hu, —u,+H;—U)] (2.18) 


8C,=coslw (4, -H—H+ Eq] 
Nestas condições [6] pode obter-se das expressões 
anteriores o resultado 


É fre : 
fui E | a [| hu DR tm, — UR HH) x 
4 
xcoslo (2: +p +Ho—Ba-H9I| | dy, + 
aí Ph, 
+— | A hapR +a-ugR tou, x 
4 
X cos[w pt E, -Ho*Hg— 44) [| du, + 
A aos 
+=] [impor +a nr tau p)x 


A — GD 


(2.19) 


4 
Xcosto (H, —-Ho- By + uq) [| du, 4 
di qa A 
+ a | ii [|] tu R tr+ Hj —ByR At+Ho— H4)X 
4 
Xcostw (2 +p,+Ho—Ho— 9] [| du, 
Para resolvermos estes integrais vamos exprimir as 


correlações R,(r) no espectro S;(w) por TF inversa, 
obtendo-se por fim 
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4 
a , J2w,l 
Ciy= 7 [altoP+ Poe D+ 


(2.20) 
ú — [BT 
+T “(ve ' +29 (U? + | JO) cos 2, 
em que se definiram 
o Es (2.21) 
Hwy)= | H(w+w JS (w)e/º! — 
— O 0 + 2n 
(2.22) 


+ qo 
a do 
H(o-wo)H (o+w)S (9) e! — 
E ai 2n 
Note-se que as 3 primeiras parcelas constituem o 
quadrado de 


dit)= | 


À — jul 
Ive O +l(ve o 
Se substituirmos (2.21) nesta expressão obtém-se 
fácilmente, tendo em conta o carácter real e par de 
[H(w)P e de Sdw) e as expressões (2.4) e (2.6), o valor 
(4/a2) RT), pelo que se tem finalmente 


d 


2 & | 9 9 (2.23) 
ClU=2R,+ B 2 | Ju) + |J(0)) cos 2w,t 


Conclui-se assim que é a presença da segunda parcela 
da expressão anterior, através da grandeza J(1) que, no 
caso geral, impede a verificação da relação (2.9), 
traduzindo assim o carácter não gaussiano do processo 
v(L). 

Notemos que a referida parcela se anulará se e só se 0 
produto Sw) H(w-wo) H(w +wp) se anular (anulando 
J(1)). Se não pusermos nenhuma restrição ao espectro 
Sw), isto equivale à condição 


H(w— w,) H (wo + w)=0 (2.24) 


Esta condição é satisfeita, dentro da classe dos filtros 
passa banda, para a situação esquemática representada 
na fig 2.3,em que w' 20 


h IH (c0)|? 


| 
| 


| 
RO 55, 


O + 0)! 3u), su) 


figura 2.3 Função de Transferência quadrática de um filtro 
satisfazendo a condição (2.24). 


Na prática tem-se a situação da fig 2.4 

Em conclusão : para filtros satisfazendo a condição 
(2.24), de acordo com a fig 2.4, verifica-se (2.9), tal como 
no caso não-heterodínico em que v(t) é gaussiano. 


2.3 A distribuição de z (t) 


A presença do factor cos(wpt + a) implica no caso geral 
que v(t) não é gaussiano, pelo que os cumulantes de z(l) 
não tomam a forma estudada no capítulo anterior, não 


o E Ff 


fo <3Fo 


figura 2.4 Situação tipica dum analisador heterodinico com filtro de 
banda estreita satisfazendo a condição (2.24). 


parecendo em princípio muito viável tirar quaisquer 
conclusões sobre a distribuição de z(t). 

Na situação porém que acabámos de estudar, em que 
é válido (2.9), o problema modifica-se. Com efeito pode 
provar-se |7l que a condição (2.24) não implica apenas 
(2.9) mas sim que ví(t) é um processo exactamente 
gaussiano. 

Nesta conformidade, não só Efz(t)J e var/z(t)/ 
apresentam os mesmos valores obtidos no cap. 1, como 
também são válidas todas as expressões (assintóticas) 
dos cumulantes aí estudados. 

Concluímos assim que, nas condições (2.24), os 
analisadores heterodínicos apresentam necessariamente 
para a saída z(t) - e correspondente estimador S/2nF) - a 
mesma distribuição que os analisadores 
não-heterodínicos, não havendo portanto distinção entre 
os 2 tipos de analisadores do ponto de vista estatístico. 
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RESUMO 


Neste artigo são analisadas as relações existentes 
entre algumas indicações espectroscópicas da energia 
electrónica na região do espaço catódico de uma 
descarga luminescente de azoto e o modelo de cinética 
electrónica num gás fracamente ionizado estudado 
por simulação Monte-Carlo. 


1 Introdução 


A análise espectroscópica de certas riscas do azoto 
possibilita a determinação de algumas propriedades do 
plasma, tais como a distribuição da energia electrónica e 
outras, permitindo considerar o seu controlo em vista de 
aplicações a laser. 

Trata-se então de um indicador útil para testar o grau 
de adequação de um modelo de cinética electrónica que é 
simulado por método Monte-Carlo e foi desenvolvido por 
J. P. Boeuf [1,2]. 

Baseado em indicações espectroscópicas fornecidas 
pelo grupo experimental do Laboratoire de Physique des 
Décharges (cf. figs. 1 e 2) o nosso trabalho consistiu em 


| (u.0.) 


BE-1 1EO 


0 2E-1 


sE-1 6E-1 
n =x/L 
figura 1 Intensidade luminosa das riscas 3914 À (em pontilhado) e 
3577 À (em continuo), ao longo da descarga luminescente 
de azoto. , 
determinar os factores relevantes na razão das 


intensidades das riscas I(3914A) / M35TTA) numa 
descarga luminescente de azoto. 

As distribuições espaciais das intensidades luminosas 
apresentam um máximo na entrada da luminescência 


Recebido para publicação em 1988. 


ABSTRACT 


This paper analyses the connections between some 
spectroscopic indications of the electron energy in the 
cathode region of a glow discharge and a model of 
electron kinetics in a weakly ionized nitrogen gas 
studied by Monte-Carlo simulation. 
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figura2 A razão das intensidades luminosas |(3914A)/1(3577Á) ao 
longo da descarga luminescente 


negativa e, como seria de esperar, a intensidade da risca 
3914A é bem mais pronunciada. 

Não obstante, a razão das duas intensidades apresenta o 
seu valor máximo antes da entrada da luminescência 
negativa e ao longo desta zona a sua razão é constante, 0 
que se revela inesperado. 


2 A Simulação pelo Método de 
Monte-Carlo 

2.1 O princípio do método 

O método de Monte-Carlo é utilizado para simular as 
trajectórias de um grande número de electrões numa 
descarga em campo eléctrico dc, incluindo os electrões 
secundários, de modo a estabelecer um balanço 
estatístico do seu comportamento e a determinar os 
valores das grandezas macroscópicas tais como as 
distribuições de densidade electrónica, velocidade de 
deriva, coeficiente de difusão, coeficiente de ionização. 
Pode-se encontrar uma exposição deste método aplicado 
à modelização da cinética electrónica na referência |2] 
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E (ev) 


Os dados iniciais da simulação Monte-Carlo são o 
campo eléctrico (fig.3), a distribuição electrónica inicial 
e as características do gás (que estão compreendidas nas 
secções eficazes de colisão electrão-molécula). 


1500 


Ç 
C 


d 


p= 0.7 Torr 
1000 + 


E (Volt /cm) 


500 | 


Oo SE-1 4 


x (em) 


figura3 Campo eléctrico entre o cátodo e o ânodo. V, = 489.5 Volt, 
d.=0.57 cm,p=0.7Torr. 


O primeiro tratamento teórico de modelização do 
espaço catódico foi feito por Von Engel e M. Steebeck [3], 
[4]. Neste modelo assume-se que a carga eléctrica perto 
do cátodo é predominantemente a dos iões que, devido à 
sua inércia, permanecem mais tempo do que os electrões 
nesta zona. A densidade dos iões n, é largamente 
superior à densidade dos electrões ne : n; > ne, 
resultando assim uma forte carga de espaço positiva e 
tal que p(x) = p, seja constante dentro da região do 
espaço catódico d,. 

À equação de Poisson 


de onde resulta o potencial 


x x 
V(x)=2V (1-5) 
cd 2d 
c c 


Naturalmente que este modelo simplifica em demasia 
a complexidade física do problema, na medida em que 
isola artificialmente o espaço catódico das outras regiões 
da descarga, em particular a luminescência negativa, de 
onde provêm parte dos iões, fotões e moléculas 
metastáveis, e que desempenham um papel essencial no 
balanço energético da descarga. Desprezamos 
igualmente as variações radiais do campo eléctrico 
axial, hipótese largamente verificada para grandes 
superfícies de descarga. Não obstante, medidas 
experimentais mostram que o campo electrico decresce 
práticamente linearmente em razoável acordo com este 
modelo teórico [5], motivo este que justifica a sua 
aplicação no presente problema. 

As secções eficazes referentes aos processos 


e+NH(x2tv=0)»e+ N! (BêL* v=0) 
2 E 2 u 


e + N,(BºIl,v=0)>. + N, (Cv =1) 
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figura4 Diagrama esquemático dos niveis de energia no azoto e 
representação de duas transições estudadas 


foram extraídas de [6] (ver figura 9). 


2.2 Resultados numéricos 
Os resultados numéricos foram obtidos por meio de 
um programa Monte-Carlo exposto em [1,2,7]. 


1) À distribuição da densidade electrónica. 

Se bem que a estatística seja baixa observa-se que os 
electrões se acumulam na zona da luminescência 
negativa, para lá da distância d. do cátodo. A frente 
súbita que surge em d, não tem significado físico mas é 
justificada pela baixa estatística. Pelo mesmo motivo, a 
densidade electrónica parece nula na zona catódica, mas 
na realidade existe um factor de multiplicação, que é o 
fluxo electrónico de, emitido pelo cátodo. 

Na zona de queda catódica, os electrões permanecem 
menos tempo porque são fortemente acelerados pelo 
campo eléctrico, (cf. fig.5). 


q 


densidade electrónica (u.a.) 


Õ | = 
O 2E-—| 4E—| 6E-1 8E-1 1EO 
n=x'L 
figura5 Distribuição espacial da densidade electrônica, em unidades 
arbitrárias 
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ii) Velocidade de deriva, V,. 
Como se pode verificar pela figura 6, a velocidade de 


velocidade de deriva (u.a.) 


6E-1 1EO 
n=x/L 
figura 6 Evolução da velocidade de deriva ao longo da descarga. 


deriva aumenta ao longo de uma zona de não-equilibrio 
devido à diminuição do campo eléctrico, mas perto do 
meio da queda catódica encontra-se quase em equilibrio, 
decrescendo posteriormente com o campo eléctrico. 

O ligeiro aumento de V, perto do ânodo deve-se à 
absorção dos electrões. 


4E-1 


BE-| 


Õ 2E-1 


iii) Energia média, e. 
Verifica-se uma grande semelhança entre a curva de 
e ede Vy, (cf. fig. 7). 


70 | 
50 


40 


e (eV) 


30 


20 + 


10 


“6E-4 


n=x/L 
figura 7 Evolução da energia média ao longo da descarga. 


O 2E—1 4E—1 1EO0 


iv) Primeiro coeficiente de Townsend, a. 

A fig. 8 ilustra a distribuição espacial de valores de a, 
o primeiro coeficiente de Townsend calculado pela 
expressão 


o, (e)vfix,e)de 


a (x)= — 


V 


X 
e normalizando a função de distribuição da seguinte 
forma: 


n, (x) = | f(x,v)dv 


onde v designa a velocidade dos electrões, e= 1/2 mv? a 
sua energia cinética e ne(x) a densidade electrónica.. 
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4E—1 6E-—1 


n=x/L 
figura 8 Evolução espacial do primeiro coeficiente de Townsend. 


BE-| 


Observa-se valores nulos enquanto os electrões não 
possuem energia suficiente para ionizar as moléculas, 
aumentando bruscamente até atingir valores 
nitidamente inferiores à distribuição de equilíbrio [2] e 
mantendo valores não nulos para lá da luminescência 
negativa, em oposição ao que seria de esperar se 
considerássemos a sua distribuição em equilibrio com o 
campo local. 

2.3 Factores determinantes na razão da 
intensidade das riscas. 


Tendo sido experimentalmente observada uma 
relação constante das intensidades luminosas ao longo 
da luminescência negativa, pretendia-se verificar por 
cálculo numérico efectuado por resolução Monte-Carlo 
da equação de Boltzmann, o grau de fiabilidade que 
oferecia esta aproximação como modelo teórico. 

A figura 9 apresenta as curvas das secções eficazes 
das transições analisadas. 


No 
| CA 
| — 2?úmm 
| y— a 
| 
| -500 0 
o 3914 
Ga 3577 E 
/ f x=08em 
| ' f f. x-12 em 
PSY 
ques ” / p= miss 
E (ev) 


figura9 Secções eficazes das riscas 3577 Á e 3914Á e função de 
distribuição electrônica em dois pontos da descarga 
luminescente 


A função de distribuição electrónica em energia, no 
início x=0.8em e no final x=1.2 em da luminescência 
negativa (cf. fig.9) mostram a existência de um número 
considerável de electrões de energia suficiente para 
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ionizar e excitar as moléculas de Ns. Esse número 
decresce à medida que se aproximam do ânodo. 

A intensidade luminosa das riscas é proporcional à 
frequência de excitação das mesmas. Estas estão 
traçadas na figura 10 a um factor 1/n. 
aproximadamente: 


+00 


no (v)=N | 


onde N é a densidade dos neutros, cálculo este efectuado 
para cada posição x da célula do espaço das fases. 


o(v)vftv)dv 


m— É) 


U.O. 


se 
x (cm) 


figura 10 Frequências de excitação e a razão (3914 Á)/1(3577 Á) no 
caso de uma distribuição uniforme (8 = 0.0). 


A razão | (3914Á) /1 (3577Á) resulta efectivamente 
constante ao longo da luminescência negativa. Foi 
utilizado neste cálculo uma secção eficaz diferencial 
simples de ionização isotrópica, tal como foi proposta por 
G.D.Alkhazov [8], e que permite avaliar a influência de 
diversos esquemas de balanço energético, tais como: um 
electrão incidente de energia cinética €,, colide com uma 
molécula e arranca um electrão de energia vu). À energia 
final pós-colisão disponível é e--ui=-V. Designemos por 
ô a percentagem de energia distribuida a um dos 
electrões. Pode-se ver representadas na fig.l1, três 
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figura 11 Primeiro coeficiente de Townsend calculado utilizando 
secções eficazes diferenciais simples de ionização nos 
seguintes casos: (1) distribuição não -uniforme; (2) e (3) 
distribuição uniforme. 


distribuições do primeiro coeficiente de Townsend a (x) 
em função da distância ao cátodo. Procurou-se 
determinar em que medida uma modificação de 6 
influiria sobre a (x) e, por consequência, sobre o factor de 
multiplicação electrónico M(x) =explJ“ a (x).dx). A 


BE—1 1EO 


curva (1) corresponde a ô = 0.2 (isto é, um dos electrões 
transporta 20% de energia disponível pós-colisão) e a 
curva (3) foi traçada para ô = 0.0 (i.e., probabilidade de 
distribuição de energia uniformemente repartida por 
tiragem de um número aleatório R equidistribuído sobre 
o intervalo [0,1] ). Note-se a igualdade de ambas as 
curvas (resultando os picos de flutuações estatísticas) e 
onde se supôs que o trabalho mínimo de extracção de um 
electrão no cátodo é em média de W, = 10 evV. 

A curva (2) foi calculada com W, = 2eV e 6 = 00€e 
apresenta uma distribuição que se mantém para além da 
entrada da luminescência negativa. O primeiro 
coeficiente de Townsend a (x) apresenta uma 
distribuição diferente daquela que se obteria para o 
coeficiente de ionização a-q(x) em campo uniforme que é 
dada por uma relação clássica do tipo 


B 


(E/p)' 
onde r é da ordem de 1/2, A e B são constantes 
dependendo da natureza do gás, E/p uma variável 
reduzida do campo eléctrico aplicado sobre a pressão do 
gás. Mas a possui valores não nulos para lá de d, devido 
à existência de um número superior de electrões de 
elevada energia em relação ao esperado para uma 
distribuição de equilíbrio, facto que é bem patente sobre 
a figura 8. 

O parametro ô não modifica significativamente a 
razão das duas intensidades (cf. fig.12). Pode-se 


1E0 PR 


A q (x) =Apexp 


| / N 
BE-1H N cul [| 
A [ Jez 
| w ! | 
BE—1 y ! | 
a RR TA 
4E—1 ) V À vp À 
PANA ANVAI 
2E—1 É dq N A 
/ | E ae VA 
e”, ( | 
e a ss 
0 2E-1 4E-1 6E-1 BE-1 1E0 
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figura 12 Variação espacial da frequência de excitação das riscas 
3914 À (curva (v,))e 3577 À (curva (vo) )e a razão (v,)(vo) 
no caso de uma secção eficaz diferencial de ionização 
não-uniforme. 
comparar no gráfico da fig.13 as frequências de excitação 
da banda B25,,' > xº is (vº=0-» v'=0) para W, = 10 
eVeW, = 2eV, verificando-se uma constrição da curva 
(2) e um acréscimo do máximo em relação à curva (1), 0 
que se deve ao aumento de probabilidade de excitação 
desta banda devido à presença de electrões mais 
energéticos. O gráfico da fig. 14 ilustra as alterações 
sofridas na razão I(3914A)/(3517Á) considerando dois 
Wo, diferentes para o cátodo, 

Da comparação da curva teórica e experimental na 
figura 15 ressalta um bom acordo, tendo em 
consideração as limitações de ordem experimental 
relacionadas com o poder de resolução da rede, a selecção 
de determinada banda espectral e sensibilidade dos 
aparelhos. 
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figura 13 Frequências de excitação directa XL, ' > B2E,* para 
duas energias de extracção no cátodo diferentes: (1) 
Wo=2eV;(2)W5=10€eV. 
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figura 14 Curva de 1(3914 À)/1(3577 À) ao longo da descarga e 
para: (1)W5=2eV; (2)W9=1Dev. 
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figura 15 Distribuições espaciais de |(3914 À)/1(3577 À). 
(i) Teórica. 
(ii) Experimental. 
Azoto/0.7 Torr/V,.=-500V/V,=0V/d,=20mm/ 
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RESUMO 

Neste artigo faz-se uma breve análise dos 
mecanismos da corrosão sulfática nas caldeiras. 

Apresenta-se uma análise termodinâmica da 
corrosão metálica em sais fundidos, em termos dos 
diagramas potencial-pSOs. Descrevem-se medições 
electrocinéticas do comportamento de aços de 
construção, para caldeiras, em meio sulfático. As 
características metalográficas da corrosão em sais 
fundidos são também analisadas, incluindo a 
carburização do substracto, os fenómenos de 
oxidação/sulferatação interna e a morfologia dos 
filmes de oxidação. 

Os efeitos benéficos de aditivos básicos para o 
controlo da corrosão são referidos e o desenvolvimento 
de um sensor electroquímico para monitorizar a 
corrosividade dos materiais em serviço é considerado. 

Aponta-se a necessidade de mais R&D para a 
minização desta forma de corrosão. 


Introduction 


High-temperature fireside corrosion in coal and oil- 
fired steam generating plants and generally in 
heat-recovery units is a well-known problem caused by 
molten inorganic sulphates deposited on superheated 
steam tubes. In modern pulverised fuel fired boilers, it is 
generally prevalent on final steam tubes (Fig. 1) which 
have fireside surface metal temperatures of 600-7000C, 


carrying steam at nominally 5680º0C/170 bar (super- 


heaters) and 5680C/40 bar (reheaters). 

The molten sulphate deposits are complex mistures of 
sodium, potassium and magnesium sulphates |4|, 
containing free sulphur trioxide which attacks the scale 
by dissolving the protective oxid as iron (and chromium) 
sulphates, for example: 


Fes03 + 38505 +3K9504=2KgFe(S0O4)s (1) 


The rate of this continuous dissolution process 
determines the rate of corrosion, which is therefore 
highly dependent upon the temperature gradient 
between the hot gas and the metal surface. Another 
important variable is the concentration of gaseous 
sulphur trioxide formed within the outer porous ash, 
since, for a given temperature, this determines the 


Original recebido para publicação em 13 de Abril de 1988. 


ABSTRACT 


Fuel ash corrosion is a well known problem in the 
Power Plant Industry, and is caused by molten alkali 
metal sulphate deposited on superheated steam tubes. 
lts consideration from a chemical and an 
electrochemical point of view Is presented. 

A thermodynamic analysis of metal corrosion in 
molten sulphates is given in terms of potential-pSOs 
diagrams. Electrokinetic measurements on the 
behaviour of constructional boiler steels in sulphate 
media are reported. Also discussed are some of the 
metalographic feactures of molten salt corrosion, 
including substrate carburization, internal oxidation/ 
sulphidation and corrosion scale morphology. 

The beneficial effects of basic additives in retarding 
this form of corrosion is outlined and the development 
of an electrochemical probe to monitor the field 
deposits corrosiveness is considered. 

It is concluded that a number ol outstanding 
questions need to be resolved [or alternative 
operational options to be considered to control this 
fireside corrosion, hopefully at an economic benefit. 


chemical stability of the molten deposit and controls the 
temperature at which dissolved metal is re-precipitated 
as solid oxide. 

The exact mechanism of the high-temperature 
fireside corrosion is a matter of controversy, but it is 
accepted that the presence of molten sulphate deposits is 
the dominant feature in the corrosion process. In 
general, it is assumed that corrosion of metals in molten 
salts is an electrochemical process [6], and therefore an 
emphasis on the electrochemical aspects of fireside 
corrosion should receive a high priority. 


Thermodynamic Analysis 


For a thermodynamic prediction of the processes 
which take place during boiler corrosion, consideration 
ofthe chemistry of four distinct phases is required; these 
phases are the metal itself, the layer of corrosion product 
overlying the metal, the sulphatic deposit formed by 
deposition ofalkali metal salts from the flue gas, and the 
flue gas itself. At temperatures below 9000€ the flue 
gases themselves contain only a small proportion of the 
S0s which would be expected on the basis of calculations 
of the equilibrium position of reaction: 


SO» + 1/2095 = SOg (2) 


This is due to the flue gases flowing too quickly 
through the boiler for equilibrium to be established. 
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Figure 1 Typical boiler, with indication of affected regions 


Some of the O, and SO are dissolved in the sulphatic 
deposit, however, and in the presence of an abundance of 
catalytically active material, 1. e. heavy metal oxides, 
equilibrium is established, the activity of the dissolved 
SO being equivalent to the calculated pressure. 


This SOg is of particular importance when 
considering the next stage of the corrosion process, 
namely the interaction of the protective corrosion 
product layer with the deposit. This may be expressed by 
the equation: 


MO +SO, =Mº* +SO;” (3) 

To a first approximation, at any given temperature it 
may be assumed that if the activity of SOg in the melt is 
greater than the equilibrium pressure of SOs which 
would be generated over the sulphate of metal M, when 
held at that temperature, then the oxide will begin to 
dissolve, and the protection it affords to the underlying 
metal will be reduced. Using thermodynamics, a further 
step can be taken in refining this concept. This is shown 
in Fig. 2, which depicts a Pourbaix diagram. Basically 
this describes the behaviour of an element, in this case 
Ni at 6000C, in contact with an 78 mole % LiçS04-13.5 
mole % K9504-8.5 mole % NagS04 melt, when exposed 
to different partial pressures of SOs and Os The enclosed 
areas indicate the thermodynamically preferred 
corrosion products for the particular sets of conditions, 
the Ni2+ region referring to NiSO4 being the preferred 
product. 

The upper log (Ni2+) scale refers to the activity of 
Ni2+ ions in solution; at zero, this indicates an activity 
of unity, so that NiSO4 is favoured, The negative 
extension of the scale indicates the activity of the Ni2+ 
ions which would exist at equilibrium, dissolved in the 
sulphatic deposit, but in a region in which NiO rather 
than NiSO4 is the preferred corrosion product. It is 
normally considered that when this activity falls below 
10-6, the metal will, to all intents, be inert. From this 
diagram it can be deduced what these conditions are for 
Ni at the given temperature. It should be appreciated 
that these equilibrium diagrams are based on simplified 
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Figure 2 Niinmolten Li,SO, + K;50, + Na;50, at 600 “C 
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assumptions and the interpretation of the metallic 
corrosion in the alkali metal sulphate melts which form 
boiler deposits, purely based on their application, may 
led to conflicting situations. In particular, the oxide 
and/or sulphide films generated at the metal-melt 
interface, do not confer a proper passivity, as in aqueous 
media [7-9]. 


Electrokinetic Measurements 


When one considers anodic dissolution of metals in 
sulphate melts, one can subdivide the metals in two 
groups: Active Metals and Active-Passive Metals, 
Materials in active condition (e.g. Ag, Cu, Ni) show a 
Tafel-like linear behaviour between the logarithm of the 
anodic dissolution current, i, and the potential, E, at the 
metal/melt interface (Fig. 3,). 

Materials which form surface films (e.g. Cr- 
containing alloys) have many similarities to Lhe active- 
passive metals in aqueous solutions (Fig. 3,b). The 
anodic polarization curves of these materials show areas 
of severe corrosion and of "passive-film" formation. The 
working temperature, the electrode potential at the 
metal/melt interface, and the acidic content of the melt 
(e.g. SOs, V9O0s5, MoO3 or other acidic compounds) 
influence the composition of the passive film, and 
therefore its ability to enhance or reduce corrosion. 

Electrokinetic measurements on a number of steels 
and other constructional boiler materials in molten 
sulphates showed that the high-chromium Cr-Ni alloys, 
of composition 23Cr-34Ni-1.25Mo-0.4Nb, 30Cr-50Ni- 
2Mo-0.2Ti-0.022Zr, 23Cr-43Ni-6W, are definitely highly 
resistant to sulphate corrosion. This satisfactory 
behaviour has also been confirmed by non- 
electrochemical observations [3]. 


Metallographic Aspects 


The major metallographic features of molten sulphate 
corrosion include carburization of the metal substrate, 
the formation of Fe-Cr-Ni spinel oxides rather than the 
protective chromium sesquioxide, and internal 
oxidation/sulphidation of the metal substrate. In our 
electrochemical studies, we did not observe any 
carburization attack, but association of the fireside 
corrosion with carburised layers is reported in detail in 
the literature [2]. 

The transformation from Cr90s type scale to spinel 
has been observed in our polarization measurements: for 
a given material which forms surface films, there is a 
critical potential above which this transformation 
occurs. This "spinel-formation" potential is a function of 
material composition, temperature, and the presence of 
oxidizing agents in the melt. In particular, SO 
stimulates the corrosion process, shifting the 
spinel-formation potential in the direction of more 
negative values and leading to the formation of 
spinel/magnetite containing substantial quantities of 
sulphides. 


C. A.C. Sequeira 
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Figure 3 - Anodic polarization curves of different materials in molten 
Li,504 + K;504 + Na;50, at 600 “€ 
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Corrosion Control and Monitoring 


Boiler corrosion can be controlled by (a) an adequate 
choice of the materials and/or coatings for use in the 
boiler, (b) the additives based on MgO or Meg(0OH)s 
directly in the boiler, (c) coal-blending, (d) low NO, 
techniques, and (e) other experiences. It is clear that 
consideration of the electrochemical aspects of this form 
of corrosion by the boiler designers, will help in the 
implementation of the measures needed for corrosin 
control. 

Boiler corrosion monitoring by electrochemical probes 
have been reported in the literature [1,5]. The 
continuous attack can be assessed by potentiometric 
electrodes which respond to the oxidizing, reducing, 
acidic and basic character of the melt. An additional 
possibility which we have been exploring is the 
assessment of a probe intended to act on boiler running 
conditions. This technique uses high electric fields and 
tries to take advantage of the conductivity of the boiler 
gaseous environment. 


Concluding remarks 


Fireside corrosion in boilers is of electrochemical 
nature. Therefore, an electrochemical approach is 
recommended. This allows the corrosion rate to be 
measured under controlled laboratory conditions during 
the course of the corrosion experiment, permiting a more 
detailed investigation of corrosion mechanisms, and this 
in turn allows the important factors which influence 
boiler corrosion to be defined, 

Our experimental data suggest that Lhe 
aggressiveness of the sulphate corrosion depends of the 
working temperature, concentration of the oxidizing 
species in the melt and the composition of the materials 
in use in boilers. For example, materials containing 
more chromium and/or less nickel show less 
sulphidation attack. 

Polarizations within the passive regions lead Lo 
different corrosion rates depending on the value of the 
polarization potential. Coupling of these polarization 
experiments with measurements of the transport 
properties of different metal oxides may lead to new 
material formulations containing alloying elements 
which improve the oxide stability. 

There are a number of outstanding questions to be 
resolved, but if a fuller understanding oÍ the overall 
electrochemical mechanism is developed this will allow 
alternative operational options to be considered to 
control this boiler corrosion, hopefully at an economic 
benefit. 
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RESUMO 

Descreve-se um método para resolução de equações 
integrais singulares em Ly(J), onde JCR é um 
intervalo, utilizando a factorização generalizada do 
correspondente símbolo. São apresentadas fórmulas 
que permitem obter essa factorização, no caso geral, e 
obtém-se explicitamente a solução de uma certa classe 
de equações integrais singulares. 


1 Introduction 


The singular integral operator of Cauchy type Su 
occurs quite often when dealing with problems in 
several areas of Physics and Engineering (such as 
Elasticity, Electrostatics, Aerodynamics) and the study 
of singular integral equations is thus of great 
importance. 

The purpose of this work is to focus attention on a 
method for studying singular integral equations which 
allows us not only to draw a conclusion as to the 
existence and uniqueness of the solution, but also to 
determine it explicitly. 

We are particularly interested in singular integral 
equations in Lg(J) of the form 


fix) 


| 
vo ->| — dx=elt) LEJ (1) 
Jt—x 


rt 
where JCR is a finite interval or the half axis 
R'=]0,+[, LJ) is the space of square-integrable 
functions on J and gE LJ) is a given function, As we 
shall see, such equations are equivalent or can be 
reduced to equations of the form 


| 4 felt x) p(x) dx = p(t) ,tERT (2) 
R 


where b,yeL2(R') and k belongs to the space of 


tempered distributionss (KR). 

If we identify in a natural way the two spaces Lo(R+) 
and Lo+(R)=(h,f, f€L9(R)) (where by h, we mean the 
characteristic function of R+), equation (2) can be 
written as 


(3) 
| ++ 
Fo = 


(4”,p" ELSR)) 
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ABSTRACT 

A method to solve singular integral equations in 
Lo(J), where JCR is an interval, using the generalized 
factorization of the corresponding symbol, is described. 
Formulas that allow us to obtain this factorization, in 
the general case, are presented and the solution of a 
certain class of singular integral equations is 
explicitly determined. 


The operator Ta, which we shall call generalized 
Wiener-Hopfoperator with pre-symbol a, is defined by 

TQo!=h Farq” , DP ELSAR) ni 
where by F we denote the Fourier transform in s' (KR), 
and a = Fk is assumed to be an element of L(R). 

Thus, to any equation of the form (1) we can associate 
a function a€ LR) (i.e., essentially bounded in R) which 
turns out to be intimately related to the Fredholm 
properties and invertibility of the corresponding 
operator. In particular, once a convenient factorization 
of the pre-symbol a is achieved, we can garantee the 
existence and uniqueness of the solution and determine 
it explicitly. 

We present here the solution of a singular integral 
equation of the form (1), with J=R+*, as well as a 
method to study this type of equations in the general 
case. 

Hinally notice that, although we have considered here 
only equations in Lgo(J), singular integral equations in 
Lo(J),I<p<e, can be studied in a similar way. 


2 Definitions and Preliminary Results 
The formula 


| 
Sab) =— | — (| tER (5) 
ni Jp x-— 


(where the Cauchy principal value is taken for the 
integral on the right-hand side, i.e. 
(x) (x) 4 
| j , dx= lim | / , dx = -—ô, t+ óI(6) 
àX— E a 
R 550 Ri e 
defines, according to a classical theorem of M. Riesz, a 
continuous operator on Lo(R). 
Since we have, for any feLo(R), 


Técnica -1/91 31 


Generalized Factorization and Singular Integral Equations in L2 


Sn belongs to the class of convolution operators in Lo(R), 
which are of the form 


T =F lar 
a 


To each of these operators we can associate a 
generalized Wiener-Hopfoperator in LR ') 


(a€L (R)) (8) 


as 


T=h T| (9) 
E “ “ILS (R) 


where a is called the pre-symbol of the operator T,. It can 
be shown (ef. [3] ) that the Fredholm properties and Lhe 
invertibility of these operators depend on the existence 
of an appropriate factorization of the corresponding 
pre-symbol, which is called generalized factorization 
relative to Lo(R). 

Denoting by Lo '(R) and Ly (R) the subspaces of Lo(R) 
consisting of the functions which are equal to zero in K 
and R+, respectively, and by L9+t and Ly the 
corresponding spaces of Fourier transforms, we shall say 
that a admits a generalized factorization relative to 

= À k 
()  a(b=a (b) e À a(o) er) 9 
= Nt+i já 
with 
REZ 
(a )JtI 
——— (E 
t-i É 
(a, ad 


É AA 
t+ - 


E =1 tio. 
(ii) (a VS la VIR ÇSLAR) 


is a bounded operator, where by Ry we denote Lhe 
subspace of Lo(R) of proper rational functions with poles 
off the real axis. 

It can be shown (ef. [3] ) that the Wiener-Hopf 
operator T, is Fredholm iff its symbol a admits a 
generalized factorization relative to Lo(R) and, in this 
case, the Fredholm index of T, is equal to -k. Morcover, 
T, is invertible iff the generalized factorization of a is 
canonical, i.e. a=aa,. The inverse operator is then 
given by 

— 1 1-1 —lm — 1 —la 
om =h F a Fh Fa F ru 

We are particulary interested in the case where Lhe 
pre-symbol is piecewise continuous and has finite limits 
at to We can define, then 


(11) 


at)+a(t”) 
Z 


a(t*)- alt” 


) 
th(my) ds 


a (tn) - 


LER, ué R) 
where define a(ot)=a( fo) and th(to)=+1. The 
function as, which is called the symbol of the operator T,, 
corresponds in the complex plane to a closed curve which 
can be obtained by adding to the image of a(t) in the 


plane the segments of straight line connecting the points 
a(t) and a(tt) whenever tis a point of discontinuity of a. 

Ifinfa pe Re as (t,H)|>0, we say that a is non-singular 
and we define ind a, as the index of the curve 
corresponding to as relative to the origin. 

In fact, the pre-symbol a admits a generalized 
factorization relative to L9(R) iff a is non-singular and, 
in this case, k = ind as (cf. [3,2] ). This provides a simple 
necessary and sulficient condition for the operator T, to 
be Fredholm (or invertible). 

lfa is non-singular, the factors a. and a 4 in (10) can 
be obtained through the identities 


(In a,) 
: é (13) 
E | É 
a (= expllt+ 1)! E pera | 

(In a,) Ina lx) 
() | 0 | 
a lt)= explt— OQ, — — | - dx | (14) 

E x—l IR I+Y 


where 


X— LO 
a )=ato( sr, 
Ú x+i 


] | 
É Ria Qr=2U-Sy) (cf 111) 

Once these factors are known we can determine the 
inverse operator (T,)-!, as indicated in (11). 

Returning now to equation (1), let us show that it can 
be reduced to a Wiener-Hopf equation in Lo(R+), 
whatever Lhe interval J may be. 

In the case where J=R+, the equivalence is obvious. 
Considering (7), it can easily be verified that the 
operator in the left-hand side of equation (1) can be 
represented as in (4), with a(t)= À-sgn t. Since its symbol 
as corresponds in the complex plane to the segment of 
straight line connecting the points z=-1| andz=+1, we 
immediately deduce that the operator is Fredholm (and 
invertible) if A€/-1,1]. The explicit solution of eg. (1) can 
be obtained as we shall see in the following section. 

As to the case where J is a finite interval, it is enough 
to consider J=/0,1]. Since 


b + 
PELOS LAR) 


— | 
Eu step Auta) LER" 


is an isometric isomorphism from L9/0,1] onto Loy(R+), 
the singular integral operator Sy defined in Lo(J) by 


l (x) 
Sfly= — | did dx (1€ J) (18) 
Hi Jg x—l 
may be identified with the operator 
am psi + 4 
M=YS,Y :LAR I>LAR) 
| (19) 
| (1/2(x— 4) 
muto =| OD qd ER 
DE JR | -explx—t) 


Now, since M can be written in the form M="T,, with 
a(t)=th(nt), tER, equation (1) is equivalent to an 
equation in Ly(R +), ofthe form 
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with a(t)=A+th(nt). 
As a final remark, notice that a similar rmethod can be 
used to reduce equations of the form 
1 | f(x) 
fo — =| desimmmeni /pE ci PRE) 
ni dy fed? 


to a Wiener-Hopf equation in Ly(R+). 


teJ(O<a<l) 


3 Example: Solution of a Singular 
Integral Equation in R+ 


Let us consider the equation 


f(x) 


l 
ro+ + | | — dx=g(b LER" (20) 
nlg” x—tl 


which, as we have seen, may be writen in the form 
DA E 

Assuming A€[-1,1], a admits a canonical generalized 
factorization a=a a, where the factors a and a, can be 


obtained from the identities (13) and (14). 
From (13) we have 


ja()=h-sgnt 2) 


Et (22) 
a ()= exp (iam) (— it) 
where 
“fa dA 
“* om É LF 


2* = expltalog 2) (2€C/H0)a €C) 


with 


log z= Injd+iarg 2 ;,—nSargz<n 

As to a., it can be obtained from (14) or, alternatively, 
from the identitya=a as: 

(23) 
a ()=(A+1)(i" 

It is easily seen that a. and a, satisfy conditions 
(1)-(Gii) in the definition of generalized factorization. Let 
us now determine the inverse operator (T,)!. 

Assuming 0O<Rea< 1/2, 


aHO=A +10" Hip? 
and 


Flaclpf=H *f o ,fELAR*) (24) 


where 


, u—l 
h M-x)(—=) (25) 


H = (141) 
4) ) Ta) 


on the other hand, 
a (0) =exp(- iam) (— i)º=exp(— iam) (io) (> in! 


where 


h lada 


(26 
(int !=FH (o) di 
(cf. [4,5]). 


From (24), (25), (26) and (11) we obtain the solution of 
equation (20): 


RH (gá—— 
e O Do) 
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fo = (Tg) 


A+1)expl= ian) d [ af 
= SS (t— x) 
PMa)Pl -a) dt 


(8 =)" !g(s)ds dx 
] x 


E sºg(s) 
| ds] (t>0) 
R + 


9 E] 

=(1—1) lAg()— — 
TIL t—s 

To obtain the last identity we had resort to a change 

in the order of integration, to the change of variable 


x=t-(t-s)/(1-y) and to the formula 


dy = cotg (an) = — ih (cf. 15]) 


q-1 
» mn —y) 
We could use a similar procedure in the case where 
l2<Rea<o. 
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RESUMO 

Este artigo descreve a concepção e o funcionamento do 
sistema CITRA - Controlo Inteligente de Tráfego para 
redes de AGVs. Dada a planta duma fábrica, a 
localização inicial de um certo número de AGVs 
("Automatic Guided Vehicles") e a localização final 
pretendida de cada um, o CITRA utiliza técnicas de 
inteligência artificial para encontrar um percurso para 
cada AGV tão curto quanto possível e que evite 
"colisões" com os outros. À necessidade de recorrer a 
métodos heurísticos deriva da explosão combinatória 
que se obtém quando se procura resolver o problema 
gerando sistematicamente todas as soluções possíveis. 


1 Introdução 


O estudo das fábricas flexíveis, ou FMSs (Flexible 
Manufacturing Systems), é hoje uma área activa de 
investigação. (Uma boa fonte de informação e de 
resultados recentes é a revista IEEE Transactions on 
Robotics and Automation). As FMSs são em geral 
compostas por um determinado número de estações de 
trabalho automatizadas, cada uma das quais executa 
uma tarefa ou um conjunto de tarefas necessária(s) à 
obtenção do(s) produto(s) final(ais). Neste contexto um 
problema importante é o do transporte de materiais, 
peças e produtos intermédios entre as várias estações. 
Uma solução é a utilização de AGVs, veículos 
automáticos cujo percurso pode ser reprogramado 
sempre que isso se justifique (p. ex.: devido a uma 
variação dos volumes de produção relativos dos vários 
produtos em fabrico, levando a uma reatribuição das 
estações; devido a modificações no processo produtivo 
e/ou no produto final, etc.). Por razões óbvias de custo e 
de eficiência interessa que o percurso de cada AGV seja o 
mais curto possível. Além disso verifica-se que, devido à 
forma como as vias de comunicação são realizadas na 
prática (p. ex. calhas ou carris), não é possível Ler mais 
que um AGV em cada troço num dado instante. É então 
necessário encontrar, dadas as estações de partida e de 
chegada, um conjunto de percursos simultâneos para os 
AGVs envolvidos que satisfaça estes dois critérios. 

O artigo é composto pela definição do problema em 
termos de entradas/saídas e de representação interna, 


Trabalho desenvolvido no laboratório da Secção de Sistemas e 
Controlo do Departamento de Engenharia Electrotécnica do IST. 
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ABSTRACT 


This paper describes CITRA, an intelligent traffic 
control system for AGV networks. Given the map of a 
factory, the inicial location of a given number of AGVs 
(Automatic Guided Vehicles) and the desired final 
location of each of them, CITRA employs artificial 
intelligence techniques to find a path for each AGV 
that is as short as possible and avoids collision with 
other AGVs. The need to resort to heuristic methods 
derives from the combinatorial explosion that occurs 
when one tries to solve the problem by systematically 
generating all possible solutions. 


que é tratada na próxima secção, pela exposição do 
método de resolução utilizado, que ocupa a seguinte, e 
por algumas considerações acerca da implementação 
realizada. 


2 Definição do Problema 


Para caracterizar completamente o problema é 
necessário definir como é representada a fábrica, como 
são representados os AGVs e os condicionalismos sobre o 
seu movimento, e a forma em que a solução é 
apresentada, 


2.1 


A planta da fábrica é especificada por meio de uma 
das listas das "ruas" existentes, sendo cada "rua" 
definida pelas coordenadas (x, y)) e (x,, y,) dos dois 
pontos P e P, que une. As ruas admitem-se horizontais 
ou verticais no plano XY, isto é, ou x =x ou y =y,, 
podendo então cada rua ser representada por uma lista 
do tipo (hyx, x,) ou (vu xy, y,) conforme for horizontal ou 
vertical respectivamente. Ás coordenadas são 
necessariamente inteiras. Por definição um cruzamento 
é um ponto P que é extremidade de 2, 3 ou 4 ruas 
simultaneamente, não devendo portanto haver ruas que 
se cortem a meio. (Por exemplo não devem existir 
simultaneamente na mesma planta as ruas (h 0-1 1) e(v 
O -1 1); essa configuração será representada pelas quatro 
ruas(h0-10),(h001),(00-10)e(v 00 1).) Para que o 
sistema seja eficiente não deve haver cruzamentos só de 
duas ruas ambas horizontais ou ambas verticais. 
Admite-se que todas as ruas podem ser percorridas em 
ambos os sentidos. 
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2.2 Movimento dos AGVs 


O ponto de partida de cada AGV é especificado pelas 
suas coordenadas x e y, e identicamente para o ponto de 
chegada. Deverá ser um extremo de qualquer rua, 
incluindo cruzamentos. Um AGV pode então ser 
representado por uma lista do tipo (x,y x, y,), em que (x, 
y,) é o ponto de partida e (x,, y,) o de chegada. Us pontos 
de partida e de chegada devem ser diferentes entre si. 
Admite-se que todos oa AGVs se movem à mesma 
velocidade, tomando-se como unidade o tempo que um 
AGYV leva a percorrer uma rua de comprimento unitário. 
O tempo total gasto por um AGV num percurso é então a 
soma dos comprimentos das ruas que percorreu com o 
número de instantes em que esteve parado. O sistema 
procura minimizar duas quantidades: a média para 
todos os AGVs do tempo gasto no percurso e o máximo 
dos tempos individuais gastos, que é o tempo que a 
operação global leva a completar-se. 

Admite-se que se dá uma colisão sempre que dois ou 
mais AGVs se encontram na mesma rua, ou seja: uma 
rua só pode ser percorrida por um AGV de cada vez, em 
qualquer sentido. Note-se que podem estar dois AGVs 
parados no mesmo cruzamento, desde que ocupem ruas 
diferentes. A situação de "troca de ruas”, em que num 
dado instante o AGV A está na rua 1 eo AGV B está na 
rua 2, e o no instante seguinte À está em 2 e Bestá em 1, 
é também considerada uma colisão. Só são consideradas 
válidas soluções em que não haja colisões de qualquer 
dos tipos. 


2.5 


O sistema CITRA procura sempre encontrar uma 
solução, mesmo que não seja a óptima. Acusará erro se 0 
problema for impossível. A solução é calculada off-line, 
sem limites ao tempo de cálculo necessário, € 
apresentada graficamente, através de animação desde o 
início até ao fim das várias operações. O CITRA fornece 
igualmente como saida utilizável por processos 
posteriores uma lista ordenada dos movimentos que cada 
AGV tem de efectuar em cada instante para ir do seu 
ponto de partida até ao ponto de chegada. Estes 
movimentos podem ser: 


Forma da solução 


eu: Para cima. 

e d: Para baixo. 

e |]: Para a esquerda. 
er: Para a direita. 

e w: Espera. 


Por exemplo, para percorrer a rua (th O O 4) seguida da 
rua (vu 4 0-3) a sequência de movimentos é rrrrddd. 

O sistema garante que nenhum dos AGVs tem 
alguma vez que recuar, isto é as listas de movimentos 
nunca conterão as sequências ulw)d, d(w)u, r(wl e 
I(w)r, onde (w) representa zero ou mais w. Apresenta 
também no final os valores do tempo médio e máximo 
correspondentes à solução adoptada. 


3 Resolução do Problema 


O primeiro passo da solução do problema consiste em 
encontrar para cada AGV o caminho óptimo desde o 
ponto de partida até ao ponto de chegada, 
independentemente dos outros AGVs. Desde que exista 
pelo menos um caminho, ele é de certeza encontrado, e se 
existir mais que um caminho possivel (que é obviamente 
o caso normal) o sistema garante que o caminho 
encontrado é de facto o óptimo, ou seja, o de menor 
comprimento . Em seguida comparam-se os percursos 
gerados para todos os AGVs, de forma a detectar se há ou 
não colisões entre AGVs, de acordo com a definição de 
colisão dada acima. Se não houver qualquer colisão, está 
imediatamente encontrada a solução do problema, e o 
CITRA devolve-a. Se houver uma ou mais colisões, O 
sistema vai então tentar resolvê-las perturbando os 
caminhos dos AGVs envolvidos nelas; ver-se-á em 
seguida em mais pormenor como isto é feito. Esta 
perturbação pode criar novas colisões, € O processo 
repete-se até ser encontrado um conjunto de caminhos 
livre de colisões. Este é então devolvido pelo CITRA. A 
solução encontrada não é necessariamente óptima, no 
sentido de minimizar os tempos máximo e/ou médio, no 
entanto, na prática tenderá a sê-lo, porque em cada 
conjunto de perturbações os caminhos individuais 
encontrados são novamente os óptimos, sujeitos as 
restrições introduzidas para resolver as colisões. 


S.l 


Para tornar mais fácil compreender como o algoritmo 
A é utilizado pelo CITRA, é conveniente começar por 
rever rapidamente o seu funcionamento. O A foi 
inicialmente proposto em [1], e mais tarde melhorado em 
[2]; a referência aqui utilizada [3] incorpora esta 
correcção. 

O A é geralmente utilizado em inteligência artificial 
para resolver problemas de travessia com custo minimo 
de "OR-graphs" = grafos em que os nós representam os 
diferentes pontos de um espaço de estados, e os ramos 
partindo de um nó representam os diferentes operadores 
que é possivel aplicar nesse ponto; há que escolher um só 
operador dentre as várias alternativas — daí o "OR". 

O A faz busca "melhor primeiro", isto é, em cada 
passo expande, dentre os nós do grafo já atingidos mas 
ainda não expandidos, o nó que parece mais prometedor. 
O objectivo é chegar à solução o mais depressa possivel, 
sem ter que fazer a expansão exaustiva de todos os nós do 
grafo. Para isso o À necessita de saber o custo de chegar 
a um nó, e de ter uma estimativa do custo de chegar ao 
objectivo a partir do nó. Prova-se [3] que se esta medida 
nunca sobrestimar o verdadeiro custo, o algoritmo 
encontra necessariamente a solução de custo minimo. 
Por outro lado, se a medida subestimar demasiado o 
custo isso levará a perda de esforço em expansões 
inúteis. 


O algoritmo A* 


| 


| Note-se que um AGYV isolado nunca Lem que parar para esperar 


que os outros passem, 
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Na totalidade o À' requer 6 parâmetros: nó inicial, 
função que calcula o custo de atingir um nó a partir 
doutro, função que calcula uma estimativa do custo 
futuro, predicado que detecta se o objectivo foi atingido, 
função que gera os sucessores de um nó, e um parâmetro 
que indica se se pretende apenas encontrar a solução, ou 
também interessa saber o caminho seguido para lá 
chegar. 


3.2 Procura de caminhos 


A procura de um caminho de um ponto inicialpara um 
ponto final através de uma determinada planta é um 
problema típico de busca em espaços de estados, em que 
o espaço de estados é o conjunto de todos os “pontos 
válidos" da planta. Em primeira análise estes são Lodos 
os pontos da planta (que são em número finito porque as 
coordenadas são inteiras e limitadas). Na prática é mais 
eficiente considerar que são apenas os cruzamentos e os 
pontos de onde os AGVs podem partir e onde podem 
chegar, ou seja: pertence ao espaço de estados um ponto 
que seja extremo de uma ou mais ruas. Isto é assim 
porque em pontos no meio de ruas só há uma alternativa 
útil, que é andar para a frente, e portanto não há 
qualquer decisão a tomar. Por conveniência, os 
operadores que fazem a transição de um ponto do espaço 
de estados para outro consideram-se como sendo as 
próprias ruas. Uma vez que não é permitido recuar, e só 
há ruas verticais e horizontais, em cada ponto há no 
máximo três alternativas, e portanto o factor de 
ramificação é convenientemente baixo. 

Dado um ponto de partida, o conjunto de todos os 
percursos possíveis a partir dele pode ser representado 
por um grafo em que cada nó representa um ponto do 
espaço de estados, e cada ramo partindo de um ponto 
representa a aplicação de um operador válido nesse 
ponto. No caso do CITRA este grafo é a própria planta, 
ou mais precisamente a parte dela que for atingível a 
partir do ponto inicial. O A é obviamente aplicável a 
esta situação. Os parâmetros utilizados são: 


e Nó inicial: Ponto de partida do AGV 

e Custo de atingir um nó a partir doutro: 
Comprimento da rua que liga directamente os dois 
pontos correspondentes. (Se essa rua não existir, a 
passagem de um nó para o outro nem sequer é 
considerada) 

e Estimativa do custo futuro; Distância que se 
percorreria se houvesse ruas ligando o nó 
directamente ao objectivo, uma horizontal para 
cobrir a diferença de coordenadas X, e uma vertical 
para cobrir a diferença de coordenadas Y. Esta 
estimativa garante a optimalidade da solução. 

e Detecção do objectivo: O objectivo foi alcançado se o 
ponto atingido for o ponto de chegada do AGV. 

e Sucessores de um nó: Pontos que estão 
directamente ligados a ele por ruas, exceptuando o 
ponto em que o AGYV esteve imediatamente antes, 
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e O nó-objectivo final é conhecido (é o ponto de 
chegada do AGV), e o que se pretende encontrar é o 
caminho para lá chegar. 

3.3 Resolução de colisões 

Achados os percursos individuais é no entanto 
necessário garantir ainda que eles não implicam colisões 
entre AGVs. Esta é de longe a parte mais complexa, pois 
envolve interacção entre os diferentes AGVs e portanto 
um conjunto de possibilidades muito maior. Encontrar 
um algoritmo de optimalidade garantida que não passe 
por uma busca exaustiva é uma tarefa extremamente 
difícil, senão impossível. O método adoptado é o 
seguinte. Considere-se um espaço de estados em cada 
elemento é um conjunto de caminhos para todos os 
AGVs, ou seja, um "candidato a solução". Podemos então 
aplicar o algoritmo A à resolução de colisões da forma 
que se descreve a seguir. 

O nó inicial é (em primeira aproximação) o conjunto 
dos caminhos óptimos gerados isoladamente para os 
vários AGVs. Pode-se definir uma medida, designada 
por "quantidade de conflito", que é calculada da forma 
que se segue. Consideremos por exemplo uma rua cuja 
ocupação desde o princípio até ao fim da "solução" foi (O 
011002222330011). A contribuição dessa rua 
para a quantidade de conflito é 2+2+2+2+3+3=14, 
ou seja, a soma para todos os instantes do número de 
AGVs que a ocupam instante após instante, excepto se 
forem apenas O ou 1 (uma vez nesse caso não há colisão 
na rua). A quantidade de conflito é a soma deste 
resultado para todas as ruas, que pode também ser 
entendida como a soma para todos os AGVs do tempo 
passado em colisão com outros AGVs. Um nó é um nó de 
sucesso (objectivo atingido) se para ele a quantidade de 
conflito for nula. Em cada passo do algoritmo interessa 
apenas expandir o nó mais prometedor, que é de todos os 
nós ainda não expandidos aquele para o qual a 
quantidade de conflito é menor. A quantidade de conflito 
é uma medida da distância do nó à solução, e 
desempenha portanto o papel de estimativa do custo 
futuro. Como interessa apenas chegar o mais 
rapidamente possível a uma solução, e a sequência de 
tentativas feita para lá chegar é irrelevante, o custo 
passado deve ser sempre nulo. 

Os sucessores de um conjunto de caminhos são os 
conjuntos de caminhos que se obtém tentando resolver as 
colisões desse caminho de todas as formas possíveis, 
Simplificadamente isto significa o seguinte. Uma colisão 
ocorre quando dois ou mais AGVs Lentam ocupar uma 
dada rua durante intervalos de Lempo que se sobrepõem 
total ou parcialmente. Se houver N AGVs envolvidos 
numa colisão podemos considerar que há N maneiras de 
a resolver, consistindo cada uma delas em reservar a 
rua para cada um dos N AGVs durante o tempo em que 
ele a pretende ocupar, obrigando portanto os restantes a 
procurar alternativas, uma das quais é de esperar até 
que a rua fique livre. (Não interessa desviar todos os 
AGVs, já que a utilização do troço em conflito faz parte 
de um caminho já óptimo, logo as alternativas são 
necessariamente de custo maior ou igual.) Para resolver 
(ou pelo menos tentar resolver) todas as M colisões de um 
dado conjunto de caminhos é necessário fazer M 
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reservas, uma para cada colisão. Como uma dada 
solução de uma colisão pode ser combinada com todas as 
soluções de todas as outras, se a colisão i envolver N, 
AGVs o número de sucessores do nó será [l,=1 . MN. O 
factor de ramificação depende portanto do número de 
colisões que se verificam, e do número de AGVs 
envolvido em cada uma delas. 

Um sucessor particular obtém-se fazendo um conjunto 
particular de M reservas, e gerando um novo conjunto de 
caminhos que tome em conta essas reservas, e que nessa 
situação seja ainda o conjunto óptimo, em termos de 
minimizar os tempos médio e total gastos. Isto pode ser 
feito pelo método já utilizado para encontrar a primeira 
tentativa, mas generalizado para: 


6 Incluir na geração de sucessores a verificação, para 
cada rua possivel, da existência ou não de reservas 
feitas para outros AGVs e que se sobreponham ao 
período de tempo durante o qual o AGV a pretende 
utilizar. 

é Incluir nos sucessores possíveis de um cruzamento 
não só as ruas que saem desse ponto, como a espera 
na rua através da qual o AGV chegou a esse 
cruzamento, até que surja um intervalo livre na 
rua seguinte, de duração igual ou superior ao 
tempo necessário a percorrê-la. Estas esperas serão 
geradas para cada uma das ruas seguintes 
possíveis de serem ocupadas. 

€ Incluir o tempo gasto em esperas no cálculo do 
custo passado, de forma a que os caminhos 
envolvendo esperas sejam tratados no A de forma 
inteiramente idêntica à dos restantes, levando à 
escolha do caminho de duração mínima qualquer 
que seja a combinação de esperas e desvios 
envolvida nele. 


Feita esta generalização verifica-se que o conjunto de 
caminhos inicial é apenas o caso particular que se obtém 
quando não há quaisquer reservas feitas. O sistema pode 
então ser simplificado considerando antes como nó 
inicial da resolução de conflitos o conjunto de caminhos 
nulos, em que portanto também não são feitas quaisquer 
reservas, e obtendo-se automaticamente a solução 
inicial como sucessor único desse nó. O funcionamento 
do sistema toma então a forma de uma chamada ao 
algoritmo A' para resolução de colisões, em que na 
função de geração de sucessores o mesmo algoritmo é 
recursivamente chamado, com parâmetros diferentes, 
para procura de caminhos. Às considerações acima 
contemplam apenas as colisões por ocupação simultânea 
de uma rua. Para evitar também as colisões por "troca de 
ruas” é necessário na detecção de sucesso verificar 
também a existência ou não desse tipo de colisões, e além 
de gerar reservas gerar também proibições, da forma 
que se segue. 

Cada "colisão por troca” pode ser resolvida de duas 
formas: proibindo o AGV A de passar para a rua 2 no 
instante seguinte, ou proibindo o AGV B de passar para 
a rua 1 nesse instante. Havendo N colisões destas haverá 
então 2º formas de as resolver simultaneamente. Cada 
uma destas alternativas tem então que ser combinada 
em cada passo do algoritmo com todas as maneiras já 
vistas de resolver as colisões do outro tipo, para gerar os 


sucessores de um dado conjunto de caminhos. Na 
geração de sucessores ao nível da procura de caminhos a 
existência de proibições deve agora ser também levada 
em conta. Seo AGV sobre o qual recai a proibição optar 
por esperar na rua em que está poderá então haver uma 
colisão por ocupação nessa rua. Esta será resolvida 
desviando dela um dos AGVs, resolvendo 
consequentemente a anterior colisão por troca. 

Resumindo, os parâmetros passados ao A' para 
resolução de colisões são: 


e Nó inicial: Conjunto de caminhos nulos (isto é, 
listas de movimentos vazias). 

e Custo de atingir um nó a partir doutro: Zero, uma 
vez que não estamos interessados na sequência 
utilizada, como já se viu. 

e Estimativa de custo futuro: Quantidade de conflito. 

e Detecção do objectivo: Quantidade de conflito nula, 

e Sucessores de um nó: Conjuntos de caminhos 
obtidos procurando resolver as colisões nesse nó de 
todas as formas possíveis introduzindo restrições 
aos percursos (no nó inicial as restrições são nulas), 

e Pretende-se apenas conhecer o nó-objectivo 
atingido, independentemente do caminho seguido 
para lá chegar. 


Obtida finalmente uma solução, ela é convertida para 
a forma de uma sequência para cada AGV de operadores 
u, d, |, re w, tal como já descrito, e apresentada 
graficamente como simulação animada. 


4 Implementação 


O sistema encontra-se implementado na linguagem 
LISP, por ser esta mais adequada que as linguagens 
convencionais à programação de algoritmos de busca 
heurística. Em particular foi utilizado o dialecto 
muLISP-86 [4], destinado a microcomputadores tipo 
IBM PC e compatíveis, o que garante uma grande 
portabilidade e facilidade de utilização do sistema. 

O sistema CITRA apresenta-se ao utilizador como 
uma função LISP, que pode ser utilizada no ambiente 
muLISP como qualquer outra função, passando-lhe os 
parâmetros adequados (planta da fábrica e AGVs na 
forma descrita, e diversas opções de funcionamento) e 
recebendo os resultados (lista de movimentos para cada 
AGV). Para o utilizar é portanto vantajoso, embora não 
imprescindível, conhecer o muLISP. Para a visualização 


da solução é necessário um PC com placa gráfica CGA ou 


superior; independentemente disso, no entanto, O 
sistema devolve sempre a lista-solução como valor da 
função. 

Informações detalhadas sobre o sistema e a sua 
utilização podem ser encontradas no relatório técnico € 
manual do utilizador [5]. 


5 Conclusão 


O sistema CITRA foi já testado num conjunto 
representativo de exemplos, tendo-se obtido resultados 
excelentes, isto é, óptimos ou sem melhoramentos 
aparentes a fazer. O tempo de cálculo necessário varia 
naturalmente com a complexidade do problema, mas nos 
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casos testados e num PC/AT é tipicamente da ordem dos 
minutos. 
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RESUMO 


Apresentam-se os fundamentos teóricos de 
cronopotenciometria e analisam-se as suas 
potencialidades para a resolução de vários 
problemas nos domínios da Quimica Moderna. São 
ainda referidos os principais factores responsáveis 
pelo desvio das relações teóricas e citadas, de modo 
sucinto, as suas mais comuns variantes. 


1 Introdução 

A cronopotenciometria é um dos poucos métodos 
electroquímicos susceptiveis de fornecer, com bastante 
precisão, resultados quantitativos acerca das reacções 
globais de eléctrodo complexas. 

O método consiste no registo, em função do tempo, do 
potencial do eléctrodo em estudo, (denominado 
usualmente eléctrodo indicador ou de trabalho), em 
relação a um eléctrodo de referência, quando da sua 
polarização por imposição duma corrente, aplicada entre 
o eléctrodo de trabalho e um electrodo auxiliar, também 
designado por contra-eléctrodo. Na Fig.l apresenta-se 
um esquema simplificado da montagem experimental 
utilizada em cronopotenciometria. 

Durante o registo do potencial em função do tempo a 
concentração da espécie electroactiva diminui (ou 
aumenta) na vizinhança imediata do eléctrodo. As leis 
da cronopotenciometria só se verificam se essa espécie 
chega à interface de eléctrodo (ou dela se afasta) 
unicamente por difusão. Esta condição implica a 
ausência Lotal de agitação e de convecção e requer a 
utilização dum eléctrolito suporte, ou indiferente, muito 
mais concentrado que a espécie electroactiva. 

Para o controle da corrente imposta ao eléctrodo de 
trabalho tem-se vindo a usar diferentes funções 
corrente-tempo [1-5] mas a técnica mais usual é a da 
cronopotenciometria a corrente constante, que é aquela 
que iremos abordar aqui. 

A forma da curva cronopotenciométrica potencial- 
tempo que se obtém através desta técnica, apresenta-se 
na Fig.2 e consiste nas seguintes regiões principais: 

Região 1 - Esta região corresponde ao estado do 

eléctrodo nos momentos iniciais da sua polarização. O 

potencial assumido pelo eléctrodo na ausência de 

polarização, ou seja o potencial de abandono, E, dá 
lugar a um potencial E,, ao qual principia a reacção 
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electroquíimica, e essa variação de potencial deve-se 
essencialmente à carga da dupla camada. 

Região 2 - Nesta região, a variação de potencial de E, 
para É, está relacionada com a reacção de eléctrodo 
propriamente dita. É, é o potencial ao qual se anula a 
concentração da espécie electroactiva junto do 
eléctrodo. 

Região 3 - Nesta região, a velocidade de descarga da 
espécie electroactiva não pode, por si só, corresponder 
à densidade de corrente imposta. Por consequência, a 
variação de potencial que se segue a É, provém da 
descarga duma outra espécie dissolvida, ou da 
decomposição electrolítica do solvente ou electrólito 
suporte. 


O tempo que decorre entre E e E, designa-se por 
tempo de transição (1). Este importante parâmetro 
introduzido por Butler e Armstrong em 1933 [6,7] 
representa o tempo ao fim do qual a concentração da 
espécie electroactiva é igual a zero junto à superficie do 
eléctrodo. 


figura 1 Esquema simplificado de uma montagem cronopotencio-- 
métrica: F - Fonte de corrente, V - Aparelho de medida e registo do 
potencial em função do tempo, 1 - Eléctrodo de trabalho, 2 - Electrodo 
de referência, 3 - Eleéctrodo auxiliar ou contra-electrodo. 

O largo espectro de aplicações da 
cronopotenciometria, abrange desde o estudo cinético 
dos processos de eléctrodo, associados ou não a reacções 
químicas e/ou fenómenos de adsorção à superficie dos 
eléctrodos, e a determinação de coeficientes de difusão de 
espécies em solução, até às diversas utilizações no 
domínio da Química Analítica, donde destacamos a 
determinação simultânea e com elevada precisão, da 
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E (Potencial do eléctrodo 
| em estudo) 


Regiao 3 
da 
Regiao 2 | 
A, 
1 Regiao 1 
E, o tempo 


figura2 Curva cronopotenciométrica tipica. 


concentração de várias substâncias presentes na mesma 
solução, sem necessidade de separação prévia, mesmo 
para teores bastantes baixos (da ordem de 10-8 M) [8]. 

De tudo isto ressalta a grande importância desta 
técnica, que tem sido, ao longo do tempo, objecto de 
diversos estudos de índole teórica por parte de vários 
investigadores, dos quais podemos salientar os de: 
Delahay [9], Reinmuth [10], Murray e Reilley [11], 
Paunovic [12], Jeager e Kuta [13], Jain e al. [5] e 
Delahay e Mamantov [8]. 

Por tal facto, no presente trabalho limitar-nos-emos a 
apresentar e discutir apenas os aspectos teóricos 
relativos aos casos de mais frequente aplicação prática. 

Nas secções que se seguem, consideramos que as duas 
condições seguintes são sistematicamente satisfeitas: 

1l- Existência de difusão linear semi-infinita, 
eléctrodo de trabalho plano e electrólito em repouso. 

2 - Existência de electrólito suporte em excesso, o qual 
permitirá o transporte da corrente, tornando os efeitos 
da migração iónica desprezáveis. 


2 Tempo de Transição 


2.1 Definição 


Considere-se que a reacção em estudo consiste apenas 
no seguinte processo simples de eléctrodo: 


O+ne 5R 
o qual é exclusivamente controlado por difusão. 

Em face das condições 1. e 2. acabadas de referir, 
verifica-se a la lei de Fick para x=0, onde x é a distância 
da espécie electroactiva ao eléctrodo, no instante t. 

Uma vez que se está a trabalhar em regime não 
estacionário, a concentração da espécie electroactiva 
para cada valor da densidade de corrente i é função do 
tempo de electrólise t e da distância da espécie ao 
eléctrodo x. 
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No caso da redução da espécie oxidante O, a relação 
entre a concentração de O, o tempo de electrólise e a 
distância ao eléctrodo é dada pela 2a lei de Fick: 


| 2 ; 
9€ lx, t) “p é Colx,t) 
at O ax? 


em que Dgy é o coeficiente de difusão da espécie O. 

A integração desta equação foi tratada por Weber 
[14], Sand [15] e Rosebrugh e Miller [16], os quais 
impondo as seguintes condições aos limites: 

a) 


Cfx,0) = C, 


no equilíbrio, a concentração de O é a mesma em 
qualquer ponto do electrólito. 
b) 


* =x tu 
lim Cotx,t) = C, 
1—+ 
a concentração de O permanece inalterável quando se 
está afastado da interface do eléctrodo. 


c) 
| fe 
Eu ) 
nk ' dx x=0 


a quantidade de espécie O consumida por unidade de 
tempo, aquando da passagem da densidade de 
corrente catódica i (que nesta exposição se 
considerará como >0) (lei de Faraday), deve difundir 
de acordo com a la lei de Fick 

chegaram à seguinte expressão para CO0(x,t) [12,17]: 


Ca a | 
C au =0"-(——)e — )+ erfe ) 
0 12 ad E 12 182. 
N 0 2D. 1 
Ú 
com 


us = | (-2)d 
ert(y)= — exp(—z') dz 
nº O 


onde z é uma variável auxiliar 

erfe(y) = 1 —-erf(y) 
Fazendo x=0 na equação anterior obtém-se a seguinte 
expressão para € (0,t): 


2;( 12 


nF(nD)!* 
através da qual é possível calcular a variação de 
concentração da espécie O à superficie do eléctrodo, com 
o tempo. 

Tendo em conta a definição de tempo de transição 
apresentada na secção 1, ter-se-á: 


CUpg=E- 


 .=0 parati =Tt 


x=0 
Donde, por substituição na equação anterior, se obtém a 
seguinte expressão para o tempo de transição 1, que não 
é mais que uma forma simplificada da equação de Sand: 
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nFC, (nD,) 
dO Ra 
21 
onde 
n  éo número de electrões postos em jogo na 
reacção electroquímica 


F  éaconstantede Faraday e 
Do éocoeficiente de difusão da espécie O. 


Esta expressão de 1!/2 evidencia que se 
representarmos graficamente 1!/2 em função de 1/i 
obter-se-á uma recta que intersecta os eixos na sua 
origem. 

Por outro lado podemos ainda dizer que ixl/2 é 
independente da densidade de corrente i, e é um valor 
constante dado pela expressão: 

i!2 = nFC? (nD,/4)“º 

Uma vez que a raiz quadrada do tempo de transição é 
proporcional à concentração da espécie O em solução, o 
método cronopotenciométrico a corrente constante, pode 
ser aplicado em Análise Química para o doseamento 
desta espécie. 

Tendo em conta o facto de que itl/? é um valor 
constante, o tempo de transição real pode ser ajustado 
por selecção adequada da densidade de corrente a 
utilizar. Por exemplo, enquanto valores muito baixos 
(inferiores a 1 milisegundos) do tempo de transição 
experimental, não são adequados, pelo facto da precisão 
da medida destes ser baixa, valores muito elevados são 
também de evitar pelo facto de virem afectados por 
eventuais vibrações na célula de electrólise, 
consequentemente, pelo fenómeno da convecção daí 
resultante. 


2.2 Determinação experimental 


Podem usar-se métodos gráficos e electrónicos para a 
determinação do tempo de transição. 

Em ambos os casos eles são empíricos em virtude da 
falta de descrições teóricas detalhadas, quer da curva de 
potencial que se segue ao tempo de transição, quer da 
curva de potencial correspondente à carga da dupla 
camada. 

É, no entanto, muito comum usar as determinações 
gráficas recomendadas por Delahay e colaboradores 
[8,18] e esquematizadas nas Figs. 3 e 4. 

O método apresentado na Fig. 3 é aplicável quando os 
troços AB e CD são quase paralelos. Mas nem sempre 
isto acontece na prática, sucedendo muitas vezes que tais 
troços são significativamente divergentes; neste caso, 0 
método exemplificado na Fig. 4 [8] é o mais adequado. 

O tempo de transição é representado pelo segmento 
UV, em que V é o ponto em que se inicia o troço linear da 
curva potencial-tempo. 

As dificuldades encontradas na determinação do 
tempo de transição pelos métodos gráficos, aumentam 
quando as ondas são distorcidas, ou mal definidas. Em 
tais casos, houve que desenvolver métodos electrónicos, 
os quais envolvem a diferenciação da onda 
cronopotenciométrica [19, 20). 
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figura3 Determinação gráfica do tempo de transição, segundo 
Delahay e Berzins [18]. 


tempo 
Determinação gráfica do tempo de transição, segundo 
Delahay e Mamantov [8]. 


figura 4 


Apresentam-se em seguida alguns estudos 
cronopotenciométricos teóricos, relativos aos modelos 
reaccionais mais comuns na prática. 


3 Reacção de Tranferência de Carga 
Reversível 


Trata-se duma reacção electroquímica do Lipo: 


O +ne =R 
Caso em que o reagente e o produto 
são solúveis 


3.1 


No caso em que a reacção electroquimica consiste 
numa transferência de carga reversível e a reacção é 
controlada por difusão, a relação potencial-concentração 
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obedece à equação de Nernst: 


Coto.) ” 


n e Es 
ColO,t) fa 
onde Eº é o potencial padrão do eléctrodo e f, ef são os 
coeficientes de actividade das espécies O e R, 
respectivamente. 
Tendo em conta a eg. de Sand, que relaciona a 
concentração de ambas as espécies à superfície do 
eléctrodo e o tempo, a equação anterior virá 


E=Eº+] 


1/2 
EP E ro 
n n 
HaDo t 
ou: 
Pato a Ei ai 
e a JR 
onde E,/4 é o potencial para t= 1/4 dado por: 
112 
DM a ds 
RT O 
ES E + ari — 
n 1 
Do 


Esta equação mostra que ÉE,;4 é independente da 
concentração da espécie electroactiva e da densidade de 
corrente, isto é: 

a? = dE, , 

alnCº dini 
Supondo que os coeficientes de difusão são iguais entre 
si, € que os coeficientes de actividade são iguais à 
unidade, tanto para a espécie oxidante, como para a 
espécie redutora, ter-se-á E,4=E'”, o que é verificado 
experimentalmente. 

Finalmente, representando E em função de (1“/1'2.1), 

obtém-se uma recta de coeficiente angular 4 igual a 
RT/nF. 


a 


3.2 Caso em que o reagente é solúvel e o 


produto é insolúvel 


Neste caso pode supôr-se que a actividade da espécie 
reduzida R é unitária, isto é CO) f.=1, então a 
equação E=f(t) pode ser recalculada, conduzindo a: 


| RT 2i RT e 
B=E | ig) == li 4 tao 
nF nF(nD,)'? *º nE 
donde: 
RT RT 12, 
E=E'+ — nec”, + — Ina - | 
nF Ú nF' Rd 
Substituindo t por 1/4 e considerando f =1, virá: 
aq BT. É 
K E +— ln — 
T/d nF 9 
donde: 
98, = E o RT 
glni alnCº  nFf 
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4 Reacção de Transferência de Carga 
Irreversível 


Agora a reacção electroquímica poder-se-á 
representar por: 


O+ne 5R 


4.1 Caso em que o reagente e o produto 


são solúveis 


A relação entre o potencial de eléctrodo e a 
concentração superficial, no caso duma transferência de 
carga totalmente irreversível, é dada pela equação 
seguinte |21]: 


Inn fK-Cº(0,t)— 


1 


= had 
an. an E 
onde a é o coeficiente de transferência de carga da 
reacção de transferência de carga irreversível, né o nº de 
electrões envolvidos na reacção global, n, é o nº de 
electrões envolvidos no passo de transferência de carga 
determinante da velocidade e KO é a constante de 
velocidade progressiva da mesma reacção. 

À equação anterior foi derivada desprezando a 
reacção regressiva, isto é, supondo que KO o» KM. 

Substituindo Co(0,t) pela expressão dada pela 
equação de Sand, e introduzindo a expressão do tempo de 
transição apresentada em 2.1, a qual é também válida 
para processos irreversíveis, Ler-se-á: 


| = 
Es RT | 2K, RT ao ao 
== - im RE + : In (1 —l ) 
an k (n Do) ! an k 
ou 
140 rh) | 
Rr MPR Re " aa 
Cane ; an FP JR 
Fa il 


O gráfico É versus In(11/2-t1/2), é à linha recta de 
coeficiente angular RT/an,l”, a partir do qual é possível 
calcular a. 

O potencial no instante t=0 depende de K 9, C0,a, ny 
e 1 pelo que o valor de K U pode calcular-se directamente 
da equação anterior desde que os restantes parâmetros 
sejam conhecidos. 

O potencial de quarto de onda E 4, será: 


() 
RT 2k, RT nFC%nD)! 
n—es 


H Se. 
anF (nD)'Z anF 41 


TA r N | 
a? RT E, RT 
———e e, = E EE 
alncº an F 
Conclui-se assim que, desde que o processo seja 
controlado pela reacção de transferência de carga 
irreversível, o gráfico E ,;q versus In Cº, ou In i, deve ser 
uma recta de coeficiente angular RT/an,F, 


dlni an h 
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4.2 Caso em que o reagente é solúvele o 


produto é insolúvel 


De acordo com a equação: 
RT 040 RT 
E = InnFK.C(0,t)-— 
an F an F 


Ini 


estabelecida por Eyring em 1939 [21], o potencial de 
eléctrodo não depende de Cp(0,t) e por conseguinte todas 
as equções mencionadas em 4.1, são também aplicáveis a 
este caso. 


5 Reacção Electroquímica em Duas 
Etapas Sucessivas 


Neste caso o modelo em estudo é: 


O, +n,e >R, 


R, +n,e >R, 

Se O, e R, se descarregam a diferentes potenciais de 
eléctrodo, observam-se duas ondas no 
cronopotenciograma. 

A primeira onda provirá da redução de O, ao 
potencial de descarga mais elevado. 

A segunda resultará da redução da espécie 
intermediária R,, mas incluirá também uma porção 
devida à redução directa da espécie O, que ainda se 
encontra a difundir para o eléctrodo. Devido a esta 
redução simultânea durante a segunda etapa, o tempo 
de transição 1, está relacionado com n e n,, pela 
seguinte expressão: 


2 O 
nF Da, Cr, | 
(— (2 n;n, + n$) 


4j 
ou então: 


= [2(ny/n) + (n/n,” | 

Equação esta que, por exemplo, foi verificada 
experimentalmente, por Berzins e Delahay [22], para a 
redução do oxigénio e do ião uranilo, e por Sequeira [23] 
para a redução do titânio em meio fundido. 

Para n=n, virá v,=3t (caso, por exemplo, da 
redução do oxigénio) [22]; para n,=2n, virá t,=8t (caso 
da redução do ião uralino)[22]: para n, = 3n, virá 
Lt,=151, (caso da redução do titânio)[23]. 

No caso em que a espécie intermediária é instável, 
isto é, em que a 1º etapa é lenta e a 24 etapa é rápida, só 
se observa uma onda, cuja duração corresponde ao nº de 
electrões global, n=n +n,. No entanto, isto não 
significa que a reacção de eléctrodo se dê numa só etapa. 


6 Duas Reações Electroquímicas 
Sucessivas 


No caso de duas reacções electroquímicas envolvendo 
substâncias diferentes, e cujas espécies oxidantes O e O, 


M.A.G. Martinse C.A.C. Sequeira 


se reduzem a potenciais diferentes: 


O, +ne 8, 


O, +ne >R, 
a primeira onda que se obtém no cronopotenciograma, 
será resultante da redução da espécie mais nobre (O ). A 
onda que se lhe segue será proveniente da redução da 
espécie mais activa (0,) e da parte da espécie O que 
ainda esteja a difundir para o eléctrodo. 

O estudo cronopotenciométrico deste modelo foi 
desenvolvido por Berzins e Delahay [22], que 
estabeleceram a seguinte relação: 


é (U2 
to, Da, 

2 12. 
(th, tg == 6 = = 
Substituindo + pelo seu valor dado pela equação de 
Sand, verifica-se que o tempo de transição 1, é diferente 
da concentração de ambos os oxidantes O e O e é dado 
pela seguinte expressão: 


E 
L pa 
2 Nai 3 9, 


Representando [ (t +, t, 2 |] em função de Co,', 
obtém-se uma linha recta, que pode ser utilizada na 
determinação quantitativa de Cg,. 

No caso em que existem q reacções consecutivas 
(q>2), e portanto q espécies oxidadas O, O, O, e q 
espécies reduzidas R,, R,, .. R ; Reilley [24] s sugeriu a 
generalização da equação anterior, para a seguinte: 


( 1/2 
n FC (1 D ) 
q qt O 
Bic TE Pa me 

1 2 q l 2 q- Mi 
para 0 caso em que todas as espécies são reduzidas a 
potenciais suficientemente diferentes, de tal modo que o 
cronopotenciograma obtido exiba q degraus. 

Neste caso, a soma dos tempos de transição de todas 
as espécies q, (T' ) encontra-se relacionada com a soma 
dos tempos de transição das (q-1) espécies, através da 
seguinte expressão: 


en,n(D. D 


N/Do Do) "Co Co, + nºD. O 


2 0, U 


os 


, 18 
o ECO (ndo ) 
Ig E 
Te TA = = 
onde 
y ae A L 


Neste caso a relação E=f(t) é dada pela seguinte 
equação[25]: 


To O=t 
RT 
E=(E) + n( 2) 
qlI2 nF (12 qi 
q=1 


onde (E ),, é 0 potencial de meia onda para o passo q. 
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O tempo correspondente ao potencial de meia-onda 
(E ),. é dado por: 


Representando (T '* E) versus Co, obter-se-á 
uma linha recta, que pode ser utilizada para a 
determinação quantitativa de Co, - 


7 Reacção Electroquímica Precedida de 
uma Reacção Química 


Neste caso, o modelo em estudo é do tipo: 
K. 
Y=0 
kt 


O+ne >R 
onde K' e K" são as constantes de velocidade dos 
processos progressivo e regressivo da reacção química 
que precede a electroquímica. 

Podemos considerar dois casos extremos para o 
tratamento deste esquema cinético. A densidades de 
corrente muito elevadas não há tempo para a ocorrência 
da transformação química e a quantidade de Y 
convertida em O durante o tempo L, é desprezável. 

A densidades de corrente baixas, o consumo lento do 
oxidante O no eléctrodo, permite a produção consecutiva 
de O, e o sistema comporta-se como se todos os reagentes 
fossem espécies O. 

Supondo que o coeficiente de difusão é o mesmo para 
ambas as espécies Y e O, o produto ir” será dado por 
[22]: 

à nECº (un 1) 12 nt2 


És. 
2 


12 
ertl K' + ke) L na 


2K(K'4 kr") 

com K (constante de equilíbrio da reacção 
química)=K'/K” ou ainda, quando o argumento da 
função erro é maior que dois, por [26]: 


nFCº (n D)!2 nU2; 


il? =—D————— 
2 12 
2x K' +) 


CO (concentração no seio da solução) será dada por: 


Do: 
C = Coftx,0) + C(x,0) 


Neste caso ir”? é uma função linear da densidade de 
corrente, podendo o valor (K'+K")” ser calculado 
através do declive desta recta. 

Um decréscimo do valor de i. “ com o aumento da 
densidade de corrente permite diagnosticar a ocorrência 
duma reacção química anterior à electroquímica. 

Para tempos de transição muito curtos, pode estimar- 
se o valor de K. Isto consegue-se expandindo a função de 
erro da equação anterior, em série. Considerando que: 


FZ 
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192 
Es kr) 2 <0.1 
resulta: 
E nFCº (mn Dj” 
[?) E ——— scans 
i=500 2(1+1/K) 


Neste caso, a quantidade da espécie O produzida por 
reacção química será desprezável, pelo que será de 
considerar apenas a respectiva concentração de 
equilibrio. 

A equação anterior pode portanto ser usada para o 
cálculo de K. 


8 Reacção Electroquímica Seguida de 
uma Reacção Química 


O modelo em estudo é, agora, do tipo: 
O+ne >R 


R “7 

Considerando que o produto £ não é reduzido nem 
oxidado na gama de potencial em que O é reduzido. 

Os critérios de diagnóstico correspondentes às 
diversas variantes deste esquema cinético, resultantes 
da combinação das reacções electroquímicas e químicas, 
com diferentes velocidades relativas, foram 
apresentados por Reinmuth [10] tendo Delahay |27] 
desenvolvido a seguinte equação para o cálculo de E em 
função do tempo, no caso em que a reacção 
electroquimica é reversível e a reacção quimica é lenta 
[5]: 


12 Km fil K, + KA 


eta folia) ori 
4 — E Dj — =— = 
2 nF 2 1+K 


1 + kk, + KZ | 


onde 1!/2 é dado pela equação de Sand. 
Para densidades de corrente muito baixas o 
argumento da função de erro é maior que 2, e portanto: 


1/2 
E =E + [to () In )| 
= + — | |n| — = 1 )=In 
12 nP (112 1+K 
Por outro lado para elevados valores da densidade de 
corrente ter-se-á: 


9 Oxidação Anódica de um Metal Seguida 
da Formação de um Complexo 


+ Ee 
M5>M!"+ne 


Mp E = Mx Porto 


Neste caso ter-se-á |27!: 
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RE (p-—n) 
iii 12 
RT MA RT (tt) 
ic Rd hi 
n | — n | 
“2a ii (1-000!2) 
X X 
Assim, representando E em função de 
un)! 
q >> 
5 : 
(1-2?) 


obtém-se uma recta de coeficiente angular RT/nF. 
Por exemplo, para p=2 e t=0, 1461, 0 3º termo 
anular-se-á e a equação anterior simplificar-se-á: 


| (p-n)- 
RT MX 
O 
nF (Co yp-lp? 
XxX. 


10 Oxidação Anódica de um Metal Seguida 
da Formação de uma Substância 
Insolúvel 


M=MY” +qe” 


Mar +qXT>MX, 


Considerando válida a aproximação simplificativa, de 
que o filme da substância insolúvel, formado no 
eléctrodo, não constitui uma barreira de, ou para O 
eléctrodo, ter-se-á, de acordo com Delahay |27]: 


qRT 2 (2 
ar ( 
nE (1/2 


onde S representa o produto de solubilidade de MX . 
Representando E em função de In[(12-4")/12] obtém- 
se uma recta de coeficiente angular qRT/nF. 


RT 
E=E +— In|S(Cº 3º 
0 nF x 


11 Processos em que Há Adsorção 


Quando há adsorção da espécie electroactiva à 
superficie do eléctrodo, a corrente através dele é 
composta de duas componentes: 

= dt 
em que iads é à corrente gerada pela descarga dos iões 
adsorvidos no eléctrodo, e a corrente ig; é determinada 
pela equação de Sand para um processo controlado por 
difusão. Nesta situação, o produto iLl/2 não será 
constante, mas crescerá com a densidade de corrente. 

Podem considerar-se dois casos limite. O primeiro 
pressupõe que o equilibrio entre a espécie adsorvida e a 
espécie em solução se estabelece antes da electrólise, e 
que a velocidade de adsorção é lenta. Neste caso, o tempo 
necessário para o estabelecimento do equilibrio 
considera-se maior que o tempo de transição 1. No 
segundo caso supõe-se que o equilibrio de adsorção é 
rápido, e que se mantém ao longo da electrólise. 
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11.1 Processo de adsorção lento 


Neste caso, o modelo em estudo é o seguinte: 


Oi E O ol 


lento 


O a +ne >5>R 


O +ne >R 
À corrente i,ds pode relacionar-se com a concentração 
superficial molar das espécies adsorvidas T, pela lei de 


Faraday: 
nFT 


Considerando que as espécies adsorvidas AR 
(adsorved reactants) e as espécies em solução SR 
(solution reactants) são espécies independentes, 
podemos considerar três modelos: 

(a) Modelo AS.SR - De acordo com este modelo, 
admite-se que os iões que primeiro se 
descarregam na interfase do eléctrodo são os que 
provêm da camada adsorvida. Após o 
empobrecimento dessa camada, verifica-se a 
electrólise das espécies em solução |[28-301. 
Aplicando este modelo, pode definir-se o tempo 
de transição por: 

= Ligo | Tai 

em que Lud: é dado pela relação acima e Lair é 
dado pela equação de Sand. A densidade de 
corrente aplicada será igual para ambos os 
processos, isto é, i=lads=Iso|. Fazendo as 
necessárias substituições, obtém-se a seguinte 
relação para o modelo AS.SR: 


L 


(nFC) nD 


41 
Esta equação mostra que representando ii em 
função de 1/1, para o modelo AR.SR., se obtém 
uma recta de coeficiente angular proporcional a 
DC? e ordenada na origem nFT. 

(b) Modelo SAR - Neste modelo supõe-se que as 
espécies adsorvida e em solução são 
electrolisadas simultâneamente durante o 
tempo Lt. À densidade de corrente aplicada será 
utilizada pelas duas descargas, e será dada por: 


m=nFET + 


E tada E oi 


As concentrações das espécies adsorvida e em 
solução anulam-se simultâneamente no fim do 
tempo 1. Portanto, pode considerar-se que: 

bad o Uol GÊ, 
obtendo-se a seguinte relação para o modelo 
SAR: 


Técnica - 1/91 47 


Aspectos Relevantes de Cronopotenciometria 


nFC (n Dy!? 
Ki=nPP+ —>——————— 


Esta equação mostra que o gráfico iL em função 
de v1/2 é linear para o modelo SAR, sendo o seu 
coeficiente angularproporcional a CD!2 e a 
ordenada na origem igual a nFT. 

(c) Modelo SR.AR - Neste modelo supõe-se que as 
espécies adsorvidas são descarregadas depois das 
espécies em solução. Por conseguinte, pode 
aplicar-se a equação de Sand à primeira etapa do 
processo (etapa controlada por difusão). 
Ultrapassado Ltajr estamos na presença dum 
processo similar ao considerado na secção 6. O 
tratamento matemático deste modelo deve-se a 
Anson [31] e Reinmuth [29,30]; 


] “sol” “ads Qi 
ii = nF "sol + dd! FF CC08 E SF 
sol us 
O tratamento aproximado deste modelo, usando 
o tempo de transição global L, deve-se a Lorenz 
[28]: 


a 


, 112 
Gu) 2 (npry2, DCD 

9;!º2 
Esta equação mostra que o gráfico (i11)!/2 versus 
(1)-!2, para o modelo SR.AR é linear, com uma 
ordenada na origem (nFI)!2 e um coeficiente 
angular proporcional a CD!/2 


11.2 Processo de adsorção rápido 


Desta vez, o modelo a estudar é o seguinte: 


O ..=0 
& si 


ad f 


rápido 


Udo ne >R 


O +ne >R 
sol 
Considerando que as duas últimas reacções ocorrem 
simultâneamente e que o equilíbrio se mantém durante 
todo o processo de redução, as características do 
cronopotenciograma dependerão da isotérmica: 

P=K (Co A 
isto é, a concentração superficial das espécies adsorvidas 
será determinada a partir da concentração da espécie O 
em solução, na interfase de eléctrodo. A expressão do 
tempo de transição foi derivada por Lorenz [28]: 


LR IR 
ae nFCtnD)* Ko mn Ko 
21 op an 
em que 
Di a 
G=ap— eric —— 
a 
K' K, 
48 


Para tempos de transição muito grandes, &P torna-se 
desprezável. Neste caso, o gráfico 1!/2 em função de 1/i, é 
linear, e a sua extrapolação para 1A=0 dá um valor 
aproximado de K, que pode ser usado para estimar & 
para pequenos tempos de transição. Se então se 
representa a expressão 


2 K p 
A picasa 
ape 
em função de 1/i, melhor se estima K, e, 
consequentemente, pode obter-se p. Este método deve 
repetir-se até obtermos um traçado linear. Um caso 
simplificado para uma isotérmica linear e para 
C(Dul2>/ é o tratado de Lorenz |28], Tatwawadi e 
Bard [32], que deduziram a seguinte equação: 


V9 4 not? nECu/D)? 
qtde) py ER 
Ra 21 
Representando C1l'2 em função de C2/i obtém-se uma 
recta de coeficiente angular proporcional a D!'2 e 
ordenada na origem igual a (n/D)!27/2, 

No estudo dos processos de adsorção, há que medir 
uma série de tempos de transição e os efeitos da adsorção 
são maximizados para tempos de transição curtos. 

Pode dizer-se, no entanto, que a cronopotenciometria 
é um método de utilização restrita no estudo de 
fenómenos de adsorção. Técnicas Lais como a 
cronocoulometria, devem ser usadas preferencialmente 
para estes estudos. Apesar disso, a cronopotenciometria 
indica-nos se um processo de adsorção faz, ou não, parte 
do processo global de eléctrodo; por outro lado, também 
fornece um valor aproximado da concentração 
superficial das espécies adsorvidas e consequentemente 
pode usar-se para indicar que alterações no sistema, Lais 
como a temperatura, a concentração, etc., podem actuar 
na adsorção. 

É ainda conveniente mencionar, que os modelos de 
adsorção descritos se referem apenas aos casos limite, 

Entre os modelos AR.SR, SAR e SR.AR, o modelo SAR 
representa o caso intermédio entre os modelos AR.SR e 
SR.AR. Por conseguinte, será mais seguro usar aquele 
modelo, sempre que o processo exibir apenas uma onda 
cronopotenciométrica, 


12 Processos em que Há Inversão de 
Corrente 


O+ne 5R 


R-ne >0 

De acordo com este modelo, a redução da espécie 
oxidante O verifica-se antes de se ter atingido o tempo de 
transição 1. Ao fim deste tempo, a corrente inverte-se, € 
a espécie redutora R, que resultará da redução de O na 
primeira etapa, oxida-se, produzindo a espécie O. 

O tratamento matemático deste modelo, (redução 
seguida de oxidação), foi desenvolvida por Berzins e 
Dlahay [22]. De acordo com a lei de Fick, e resolvendo o 
problema das condições aos limites, em termos de Cp(x, 
tI) e Co(x, t!), onde t] é o tempo correspondente à 
electrólise inversa (oxidação), foi obtida a seguinte 
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expressão para tl, para 0 caso em que a espécie R é 
solúvel [5,22]: 


2 
d=|—— | 
| 2 
(0+1!) — 9? 
o: 
= — E E 
nb, nED 


onde os termos 1 e i se referem à etapa electroquímica 
progressiva (redução), tl e il dizem respeito ao processo 
electroquímico inverso ou regressivo (oxidação) e Dk 
designa o coeficiente de difusão da espécie R. 

Quando i e il são iguais, as funções O e À! são também 
iguais, podendo escrever-se: 


Esta relação pode usar-se como um critério de 
diagnóstico da solubilidade da espécie R no electrólito. 

Se R é insolúvel e permanece no eléctrodo, obter-se-á 
[10]: 


= 
Quanto às expressões que traduzem o potencial em 
função do tempo [5], dependem do grau de 
reversibilidade do processo de eléctrodo, sendo diferentes 
conforme este é irreversível, reversível, ou 
quase-reversível. 


13 Outros Estudos Cronopotenciométricos 


Para além dos tratamentos teóricos apresentados, 
outros esquemas reaccionais foram estudados pela 
técnica da cronopotenciometria a corrente constante. 
Dentre estes e sem pretendermos ser exaustivos, 
podemos destacar os seguintes estudos: 

-Reacção electroquimica seguida por um processo 

catalítico [47], o que pode ser esquematizado do 

seguinte modo: 
O+ne >R 
R+27z >0 

-Reacções de redução consecutivas, seguidas por 

inversão, ou cessação da corrente [33, 34]. 

“Reacção de oxidação seguida por reacção de redução, 

em que a espécie reduzida é solúvel e se encontra 

envolvida numa subsequente reacção química de 1º 
ordem [35|. 

-Idem em que a espécie reduzida é insolúvel [36]. 

-Reacção de oxidação seguida de redução, em que a 

espécie reduzida sofre decomposição imediata e 

irreversível, em duas ou mais espécies, as quais não 

são electroquímicamente activas ao potencial da 
reação directa, mas reagem após inversão da corrente 

[5, 371. 


14 Desvios das Relações Teóricas 


O desenvolvimento das equações cronopotenciomé- 
tricas, para qualquer dos possíveis esquemas de reacção, 
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é efectuado Lendo em atenção as condições postuladas 
para os sistemas ideais. 

Experimentalmente, é necessário fazer as medições 
na zona em que estas condições são aplicáveis e evitar 
interferências nos fenómenos que interessa estudar, O 


que nem sempre se torna possível. E assim importante 


conhecer os principais factores responsáveis pelos 


desvios entre os resultados experimentais e os 


correspondentes teóricos, pelo que os iremos analisar de 
seguida. 

14.1 Não linearidade do campo de 
difusão 

Deve usar-se sempre um eléctrodo plano para evitar a 
não linearidade do campo de difusão e as condições 
experimentais devem ser tais, que a difusão seja 
uni-direccional e perpendicular à superficie do eléctrodo. 

Quando isso é impossível, passa a verificar-se a 
difusão cilíndrica, ou esférica, (caso, por exemplo, dos 
eléctrodos de metais líquidos). 

Estas condições foram estudadas por diversos 
investigadores [5, 38-40], que sugeriram os seguintes 
parâmetros correctivos para a equação de Sand: 

“Para a difusão cilindrica 


E CAD)? 
i 2 — | 


º (nD 9? | 


-Para a difusão esférica 


1/2 


nFC%(nD )'? | 
Ro dd RA: 


| | 
2 (nb a 
0 
| — ——————— 
2r 
Nestas relações r é o raio do cilindro ou da esfera. 

Se as condições usadas são tais, que o termo 
(nDr)!'2/zr é muito pequeno relativamente à unidade, ou, 
por outras palavras, se as dimensões do eléctrodo (r) são 
grandes comparadas com a espessura da camada de 
difusão, dada por (Di)!2, então a equação de Sand é 
directamente aplicável sem qualquer correcção. 


14.2 Migração 


Na ausência de um electrólito suporte, o tempo de 
transição pode ser incrementado pelos efeitos da 
migração, o que foi primeiramente observado por Sand 
[15], em 1901. 

Nessas condições a densidade de corrente de descarga 
virá multiplicada por um termo correctivo (1-L'), onde t' 
representa o nº de transporte catiónico. Neste caso a 
equação de Sand virá: 


nE CD )'2 
po ir 

2(1-t') 
Na prática esta equação não é inteiramente válida, pelo 
que estudos posteriores |41| conduziram ao 
estabelecimento de expressões mais rigorosas, e sua 
verificação experimental. 
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14.3 Convecção 


A convecção consiste num movimento da solução em 
relação ao eléctrodo, que pode ter diferentes origens, 
como por exemplo a transmissão acidental de vibrações à 
célula de medida, ou a existência dum gradiente de 
densidade provocado por uma eventual diferença de 
densidades entre as substâncias que se produzem e as 
que se consomem no eléctrodo. 

Este fenómeno é particularmente significativo no caso 
dos eléctrodos colocados em posição vertical [42]. 
Orientando os eléctrodos horizontalmente, de modo que 
as espécies menos densas permaneçam sobre as mais 
densas, obtém-se uma minimização de tal fenómeno. 

Pode ainda diminuir-se os efeitos da convecção pela 
utilização de eléctrodos munidos de uma protecção [17]. 

Esta protecção apresenta a vantagem de orientar os 
gradientes de concentração segundo uma direcção 
normal à superfície do eléctrodo, contribuindo deste 
modo para a aproximação das condições ideais de difusão 
linear. Utilizando eléctrodos protegidos e com uma 
orientação adequada, pode-se manter as condições de 
difusão linear, durante cerca de 5 min. 

Os efeitos de convecção são mais evidentes para 
tempos de transição grandes e resultarão num 
incremento do produto ic!/2, pelo que convirá usar 
tempos de transição curtos para assegurar a 
insignificância de tais efeitos [43, 44]. 


14.4 Efeito da dupla camada 


Quando se aplica a um eléctrodo um impulso de 
corrente constante, uma fracção dessa corrente é usada 
para carregar o condensador constituido pela dupla 
camada electroquímica ao potencial da reacção de 
eléctrodo e só a restante é utilizada para o desenrolar de 
tal reacção, isto é: 

dé b 
onde i, é a corrente capacitiva e ira corrente farádica. 

O efeito da carga da dupla camada traduz-se na 
prática na distorção da parte inicial da curva 
potencial-tempo. À avaliação de tal efeito foi efectuada 
por Gierst e Juliard [45] e Bard [42], supondo que a carga 
da dupla camada é constante durante o tempo 1. Ter-se-á 
então: 


q, = C. AE 
onde C, é a capacidade da dupla camada, donde: 
AE 
L=C — 
ho E L 
Então: 
 nFC DP AE 
CC aue Crer 
po MNFCmD ap 
Lt = ———e. pf), ce 
9 e JR 


Esta equação mostra que para uma dada concentração, a 
influência da carga da dupla camada será mais 


pronunciada para tempos de transição curtos, em que o 
termo "C.AE/1'”" será maior, que para tempos de 
transição grandes. 

Para tempos de transição muito grandes pode mesmo 
dizer-se que a carga da dupla camada não provoca 
qualquer desvio. 

Considerando que toda a corrente aplicada, i, está 
envolvida na etapa electroquímica, ter-se-á: 


Q,=1 L 
e 
Q., “0%, AE 
IL 
€ 
ou: 


q, 4iC SE 


“, np? (Cc) nD 
onde Q. é a quantidade de electricidade usada para 
carregar a dupla camada e Q, a quantidade de 
electricidade utilizada na reacção de eléctrodo. 

Esta última equação mostra que a minimização da 
razão QQ,, isto é, a redução dos efeitos da dupla 
camada, implica o uso de concentrações (C") elevadas e 
baixas densidades de corrente (i). 


14.5 Rugosidade superficial do eléctrodo 


No caso de elécirodos rugosos, a difusão não é 
estritamente linear e a equação de Sand não pode ser 
aplicada nessas condições. 

Considerando as irregularidades superficiais do 
eléctrodo como um acidente de perfil sinusoidal, haverá 
uma contribuição esférica para a difusão. 

Supondo que o eléctrodo funciona em regime 
estacionário, Wagner [46] mostrou que, para uma 
espessura da dupla camada grande relativamente à 
amplitude da rugosidade, a área de difusão é idêntica à 
área geométrica do eléctrodo. 

À medida que a espessura da dupla camada decresce, 
a rugosidade do eléctrodo torna-se mais significativa, e a 
área de difusão aproxima-se da área superficial real do 
eléctrodo. São então de esperar desvios positivos do 
produto it!/2, à medida que a densidade de corrente 
aumenta, em virtude da rugosidade do eléctrodo. 

O efeito da rugosidade superficial do eléctrodo, foi 
ainda discutido por vários outros investigadores, dos 
quais salientamos Gerischer e Delahay [47], Miller e 
Hackerman [48], Reinmuth [29] e Bard [42]. Este último 
sugeriu mesmo, em 1963 [49], um método para corrigir 
os efeitos da rugosidade do eléctrodo, assim como os da 
adsorção à superficie deste, oxidação e carga da dupla 
camada, nas determinações do tempo de transição. 


14.6 Adsorção de substâncias não 
electroactivas à superfície do 
eléctrodo 


Em geral, a adsorção de impurezas causa um 
bloqueamento parcial da superficie do eléctrodo 
resultando num aumento da densidade de corrente e, 
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consequentemente, num decréscimo do tempo de 
transição L. 

Gierst e Juliard [45] observaram um decréscimo 
linear do produto it!/2? com a densidade de corrente ao 
adicionarem gelatina a uma solução contendo cadmiões. 

Outros efeitos, tais como uma alteração dos 
párametros cinéticos, (por exemplo,decréscimo na 
constante de velocidade e no coeficiente de 
transferência) e a introducção de uma nova etapa lenta 
no processo electroquímico (transferência de reagentes 
através da camada adsorvida à superfície do eléctrodo), 
foram também mencionados por Paunovic [12]. 


14.7 Películas de óxido na superfície do 
eléctrodo 


O efeito nas curvas potenciométricas, da formação ou 
redução de um filme de óxido, é semelhante ao efeito da 
adsorção de espécies electroactivas à superfície do 
eléctrodo: aumento de itl/2 com o aumento de densidade 
de corrente. 

A formação de tais películas pode preceder, seguir, ou 
ocorrer simultâneamente, com a reacção de eléctrodo. 

Anson e Lingane [50] e Bard [42, 51], estabeleceram 
um factor correctivo aser aplicado neste caso na equação 
de Sand, factor esse que contém a quantidade de 
electricidade envolvida no processo (Q). Após esta 
correcção, a equação de Sand apresentará o seguinte 
aspecto: 


nFC (nD)!2 
9 


pn 2. 
(142 


15 Aplicações da Cronopotenciometria 


Os principais domínios de aplicação desta técnica são 
os seguintes: 


15.1 Estudo do mecanismo de reacções 
electroquímicas e cinética dos 
processos de eléctrodo 


Reinmuth, em 1960 [10], propôs quatro critérios de 
diagnóstico, designadamente: 

1) O primeiro baseia-se no traçado de uma 

representação linear de uma função logaritmica do 

tempo t e do tempo de transição x, In G (t, x) versus o 

potencial E 


E= fllnG (tr)] 
e no cálculo do declive desta recta, representado 
matemáticamente por: 


9f[In G (tv) ] 
aE 

ul/Z (1/2: 1/2. 41/2/ç1/2, (12 - 41/2)/,1/2 são exemplos da 
função G (t, v). 

2) O segundo baseia-se na razão dos tempos de 
transição da reacção directa e inversa, determinados 
pelo método cronopotenciométrico com inversão de 
corrente: 


L/L 
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3) O terceiro baseia-se na variação de E,4 com o 
logaritmo da densidade de corrente imposta: 
dE 
dln i 
4) O último baseia-se na variação de E,4 com o 
logaritmo da concentração no seio da solução, da 
espécie que reage à superficie do eléctrodo, Co: 
E, 
aln Cº 
Aplicando simultâneamente estes quatro critérios, é 
possível determinar, pelo menos nalguns casos, qual o 
verdadeiro mecanismo da reacção em estudo, o que 
iremos a seguir exemplificar:- 
Assim, se considerarmos G (t,t) = 1/2 - t1/2/71/2, 
representarmos o potencial E em função do InG(t,t) e 


obtivermos uma recta de coeficiente angular igual a 


RT/nF e ainda: 
. O. 
dE !9inC = 0 
96 /ôtni = O 


v/v = 1/3 
então podemos concluir que se trata de uma reacção 
electroquimica reversível, que pode ser representada 
esquemáticamente por: 


O+ne EK 
onde O e R designam respectivamente as espécies 
oxidada e reduzida as quais neste caso são ambas 
solúveis. 

Contudo, para além destes, outros critérios de 
diagnóstico se podem extrair do exposto anteriormente, 
para a determinação adequada dos diversos esquemas 
cinéticos. De entre vários trabalhos publicados sobre 
este assunto [50,52,53] podemos destacar um de 
Reinmuth publicado em 1960 [10]. 


15.2 Determinação de coeficientes de 
difusão 
Tendo em conta a expressão simplificada da lei de 
Sand, que permite calcular 11/2 designadamente: 


1/2 
” nFCo (nDo) 
paia 
21 
e representando it!/2 em função de i, obter-se-á uma 
recta paralela ao eixo das abcissas, uma vez que itl/2 é 
independente de i, e que intersecta o eixo das 
coordenadas num ponto dado por nFCo“(nDç)!/2/2 donde 
é possível calcular Dq tendo em conta que nFCç e n são 
valores constantes e conhecidos. 
Um segundo processo susceptível de ser utilizado para 
o cálculo do coeficiente de difusão Do é através do 
coeficiente angular da recta de irl/2 em função de C, o 
qual é igual a 


nF(nD,)'? 
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Senderoff e Mellors [52,54] fizeram ainda 
determinações do coeficiente de difusão Dy a diferentes 
temperaturas, e através da representação desses valores 
em função de i/T, caleularam a energia de activação para 
a difusão das espécies reagentes. 


15.3 Estudo dos iões complexos 


O estudo do mecanismo de redução de i0es complexos 
e a cinética de formação e dissociação desses 10€5 
complexos pode ser, por exemplo, efectuada através dos 
conceitos teóricos desenvolvidos nas secções 5 e 10, 

Tais estudos têm vindo a ser efectuados por vários 
investigadores, de que se destacam Delahay,|18], Inman 
e Bockris [55]. 


15.4 


Laitinem e outros [56] estudaram a aplicação do 
método cronopotenciométrico ao estudo quantitativo da 
adsorção de espécies a superficie do elétrodo e 
concluiram que à presença do fenômeno de adsorção na 
gama de valores de potencial correspondente aos tempos 
de transição medidos, provoca uma dilatação dos tempos 
de transição, assim como um cocficiente de difusão 
dependente da concentração e de valor mais elevado que 
na ausência deste fenómeno. Verifica-se ainda que 
adsorção de espécies electroactivas à superficie do 
electrodo provoca um aumento de i.!/2 em relação ao 
teórico e o acréscimo deste valor com o aumento da 
densidade de corrente. 

O cálculo quantitativo da concentração à superficie do 
electrodo, da espécie adsorvida I requer o conhecimento 
prévio do mecanismo da reação do electrodo, uma vez 
que a equação i=f(1) através da qual é possivel calcular 
[', varia com o mecanismo da reacção. 


Estudos de adsorção 


15.5 Doseamento de espécies em solução 


(análise química) 


Representando ivl/2 em função da concentração da 
espécie electroactiva e desde que o respectivo coeficiente 
de difusão não varie, obtem-se uma linha recta que pode 
ser usada como uma recta de calibração, para a 
determinação do valor de tal concentração em solução. À 
precisão da determinação da concentração por este 
método, depende da ausência de variação de corrente e 
da constância de temperatura durante o ensaio, asssim 
como da precisão da medida do tempo de trasição. 

Uma vez que a linearidade do traçado de 1112 em 
função da concentração, exige a constância do 
coeficiente de difusão, a flutuação de temperatura que 
pode ser tolerada durante os ensaios é determinada pela 
taxa de variação do coeficiente de difusão com a tempera 
tura. 


16 Outras Técnicas 
Cronopotenciométricas 


Para além da técnica aqui analisada, que é a classica 
cronopotenciometria a corrente constante, existem 
outros métodos cronopotenciometrico dos quais podemos 
salientar os seguintes; 


- Cronopotenciometria com corrente programada, em 
que a corrente imposta ao electrodo varia em função 
do tempo, linearmente, exponencialmente, com uma 
potência, raiz quadrada, ou raiz de ordem n, ou ainda, 
como uma função sinusoidal do tempo, 

- Cronopotenciometria de corrente alternada (em que 

a corrente continua de electrólise se sobrepõe uma 

corrente alternada, normalmente de pequena e 

constante amplitude). 

- Cronopotenciometria derivada (baseada na medida 

da velocidade de variação de potêncial É com 0 tempo, 

representada pela |º derivada dk/dt, em função do 
tempo). 

Embora seja uma técnica eleborada, a 
cronopotenciometria com corrente programada é 
particularmente vantajosa para o estudo de sistemas de 
vários componentes e reações plurielectrónicas [17]. 

Quando a cronopotenciometria derivada, tal como a 
cronopotenciometria de corrente alternada, permite uma 
medida de tempo de transição com bastante maior 
precisão, que a cronopotenciometria a corrente constante 
[5]. 

Contudo, apesar das suas vantagens incontestáveis, a 
cronopotenciometria derivada apresenta alguns 
inconvenientes. Por exemplo, para a aplicação deste 
método, é necessário possuir informação prévia sobre a 
reversibilidade da reacção electroquimica do processo 
em estudo e no caso de tal reacção ser irreversível, 
conhecer o valor do coeficiente a. 

Finalmente, podemos ainda citar a 
conopotenciometria cíclica, como uma variante da já 
referida cronopotenciometria com inversão do sentido da 
corrente, e na qual a corrente aplicada é sucessivamente 
invertida, após cada tempo de Lransição. 
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RESUMO ABSTRACT 


Usando as relações directas entre os coeficientes do Using direct relations between linear predictor 
preditor linear e o cepstrum, descobertos por Schroeder coefficients and cepstrum discovered by Schroeder [1], 
[1], um algoritmo não-recursivo para factorização a nonrecursive algorithm for spectral factorization is 
espectral é proposto. proposed. 


1 Introduction os 

a(2!)= 2ar a aa 
In a later paper [1] Schroeder showed how to relate is dim 

directly cepstrum and predictor coefficients. We are 


; m =, 
going to describe those relations. ( p=” ] o E (8) 
Let A(z) be a minimum phase polynomial a Ri EE 2 C,2 E 2 Ca 


N This fact will be useful in the next section. 


india = RI (1) 
(2) Esp “q 2 Cepstral factorization 


and C(z) the corresponding cepstral Z-transform: 


Let R(z) be given by 
co 9) 
j (2) E SE ( 
C(z) = > Ga R(z)=A(z). A(z ) 
| 
Schroder showed that the cepstrum can be computed cd E (10) 
from es 4 rz 
-N 
(a thgthgt th 1) 7 k, | ky (3) with 
en 0) | am) 
IT N N-n 
with n= 1,2,...and the sum is to be taken over k (i=1,..., qa am (11) 
N) such that n i itn 


N In order to apply the previous results we have to 
>. E Ef (4) —transform R(z) into a minimum or a maximum phase 
' polynomial. Rewrite (10) as: 


i=1 


On the other hand, predictor coefficients, can be N 
mputed from the cepstrum by: | (12) 
compu R(z) = =. ME — 4a") 
E 1 k, / ad 0 
a = > ri ÃO ) À rs. ) i (5) Let T be a lower triangular matrix such that: 
é À pe ! l n 
a | 
with n=1,..., N andthe are such that e 
I z+z 
n (6) z+27! E: 2 
3 lk,=n | 2427? [2+z 
i=1 Co ler) (13) 
It is not hard to conclude that the cepstrum 
polynomial corresponding to a maximum phase 
polinomial A (z'), EA ad ay 
Z+2 


This matrix is easily obtained recursively and its 
Original entregue para publicação em 6 de Dezembro de 1990 elements are given by: 
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For example with n=6 we obtain: 


0 I<] 
TAF ssirts iz), i>2,]>0 (14) A É 
ij Ea lá o 01 
=D.20 j=0,1=2 2 0 1 
with To=T,=landT = O, | RR 
For example, ifn=6, we have: O 401]. 
1 O 0 O IO Ds tro 
O 1 20 0 150601 
20 1 with this transformation it is a simple task to get C'(z) 
| | | from CY(w): 
=4 0 [. | Às 
2 0 -4 0 TI 
O 5 0 -5 010 =p 
=2 0 9 0 =BD1I 
Using (14) and with the change of variable defined by 
15 
w=2+2"! 15) (19) 
R(z) is converted in Feu 
FL 
E 16 
R(w)= > yu" US) 
L=Ú 
were 
we have: 
LT é” E] Fo qo 
CHa)=> elf =5 leva "= S crer (20) 
T n “E n a n" 
= 7. (17) The use of (5) allows us to obtain a fn= 1... N). Itis 
not hard to conclude that the algorithm remains valid if 
all the zeros of A(z) are on the unit circle, but are not 
Sm Ea real. 


The algorithm we are going to describe is valid only if 
the system is minimum or maximum phase. We have to 
find the root distribution. If R(w) is minimum or 
maximum phase, we can apply it. 

We do not lose generality in normalizing y, to 1. This 
leads to the possibility of a direct application of 
Schroder's formula (3) to obtain the cepstrum C?(z). To 
obtain C'(z), the cepstrum corresponding to R(z) we have 
to use the inverse T“ of T. T is recursively computed 
and it is given by: 


0 i<j 
= NA | Ed | si (18) 
EG E) Cesijot  Tentgas 1520 
=] E Ri go 
Ma j=0,i>0 
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3 Conclusions 


The algorithm we have just proposed may be an 
alternative to Lhe usual factorization algorithms at least 
for values of N not very high. However, the 
computational cost may be too extensive in high order 
values 
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